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ABSTRACT 

Aims. We aim for a better understanding of the imprints that the nuclear activity in galaxies leaves in the molecular gas. 

Methods. We used the IRAM 30m telescope to observe the frequency range ~[86-116] GHz towards the central regions of the star- 
burst galaxies M83, M82, and NGC253, the galaxies hosting an active galactic nucleus (AGN) M51, NGC 1068, and NGC7469, 
and the ultra-luminous infrared galaxies (ULIRGs) Arp220 and Mrk231. Assuming local thermodynamic equilibrium (LTE), we cal¬ 
culated the column densities of 27 molecules and 10 isotopologues (or their upper limits in case of non-detections). 

Results. Among others, we report the first tentative detections of CH 3 CHO, HNCO, and NS in M82 and, for the first time in the 
extragalactic medium, HC 5 N in NGC 253. Her recombination lines were only found in M 82 and NGC 253. Vibrationally excited lines 
of HC 3 N were only detected in Arp 220. CH 3 CCH emission is only seen in the starburst-dominated galaxies. By comparison of the 
fractional abundances among the galaxies, we looked for the molecules that are best suited to characterise the chemistry of each group 
of galaxies (starbursts, AGNs and ULIRGs), as well as the differences among galaxies within the same group. 

Conclusions. To characterise and compare starburst galaxies, suitable species are CH 3 OH and HNCO as tracers of large-scale shocks, 
dominating early to intermediate starbursts stages, and CH 3 CCH, c-C 3 H 2 and HCO as tracers of UV fields, which control the inter¬ 
mediate to old/post starbursts phases. M 83 shows signs of a shock-dominated environment. NGC 253 is characterised by both strong 
shocks and some UV fields. M 82 stands out by bright PDR tracers, pointing to an UV field-dominated environment. Regarding 
AGNs, the abundances of HCN and CN (previously claimed as enhanced in AGNs) in M 51 are similar to those in starburst galaxies, 
while the HCN/HCO + ratio is high in M 51 and NGC 1068, but not in NGC 7469. We neither find a correlation between the HCN/CS 
ratio (recently claimed as a possible starburst/AGN discriminator) and the AGN activity. However, high enough spatial resolution to 
separate their circumnuclear disks (CNDs) from the surrounding star-forming regions is needed to find molecular abundance trends 
in AGNs. High abundances of H 13 CN and HC 3 N, as well as a similarity between the column densities of 13 CO and C 18 0, are repre¬ 
sentative of the molecular interstellar medium (ISM) in the ULIRGs. Furthermore, the chemistry of Arp 220 points towards a more 
starburst-dominated environment, while that of Mrk231 resembles more the AGNs of our sample. 

Key words. ISM: molecules - galaxies:ISM- galaxies: nuclei - galaxies: active- galaxies: abundances- radio lines: ISM 


1. Introduction 


The creation and destruction of molecular species in the dense 
ISM are directly related to a number of physical mechanisms, 
namely large-scale shocks, cosmic rays, UV radiation fields and 
X-rays. Hence, studying molecular abundances can be a way to 
infer the importance of these processes for the star formation, 
and to understand how they favour dispersion or gravitational 
collapse of compact clouds or their ionisation. In the particular 
case of the dense gas found to surround galactic nuclei, one can 
even attempt to deduce the nature of their powering sources us¬ 
ing millimetre(mm)/sub-mm spectra (e.g. lMartln et al.f201 ill . 

Comparing chemical models with observations, it has been 
claimed that the abundances of some species, such as HCN 
and CN, are enhanced in the presence of X-rays and cosmic 
ray ionisation fields that are dominant in AGNs dKohno et al 


2001 


2013 


Meiierink et al] 120071 : iKrips et ahl 120081: lAladro et al 
llzumi et alJl2013h . On the other hand, less energetic radi¬ 


ation fro m UV fields destroys heavy molecu les like HNCO and 
CH 3 OH (Mar tin et alJl2006al I2008L l20 09a). while at the same 


time helps creating others that are formed through ion-molecule 


reactions, such as c-C 3 Ht, CO + or HOC + (e.g. 

Savage,& Ziurys) 

12004; Martin et aD[2009a; Aladro et aDl20T lb 

). Shocks caused 


by stellar winds, supernovae, or molecular cloud collapse can 
release species formed on dust grains to the gaseo us phase 
(e.g. CH 3 OH , HNC O, or sulfur-bearing specie s, mti et al.ll2001b 
iMartln et al.ll2009at iJimenez-Serra et al.ll2008h . 


Unbia sed extragalactic m o lecular line surveys at mm wave¬ 
lengths (MartmetaL 2006b _ MuUeretaL^ 20L1: l Aladro et al.l 
12011 hi 1201 3l iDavis et all 1201 3l IWatanahe et all 1201 4) allow 
us to observe many species simultaneously and identify those 
that provide the best information about the physical properties 
around the nuclear regions of active galaxies. Some molecules 
much less abundant than the commonly observed species, CO, 
HCN or HCO + , seem to be more sensitive to their environ- 
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ments, and show larger differences in their abundances among 
galaxy types. They can be even used to identify the evolution¬ 
ary stage of galaxies. For example, the HNCO abundance varies 
by nearly two orders of magnitude between young and old star- 
bursts, and thus appears to be an excellent diagnostic tool of their 
evolutionary state (iMartln et al.ll2009al) . CH 3 CCH, on the other 
hand, is one of the most abundant species in starburst galaxies 
(even more than HCN and HC O + ), but has never been de tected 
in any AGN-dominated source (lAladro et al.ll201 lall2013l) . Yet, 
extragalactic line surveys carried out so far have focused exclu¬ 
sively on_a_single_type_ofgalaxy, either starb ursts (Martin etalj 
lUOTTbf). (U)LIRG s (ICostadiola et al.ll 201 It 
Davis et aD l2013h . or AGNs (I Aladro et al] 


2006b; Aladr o et al. 

Martin et al. 2011 


20131) . A detailed and uniform comparison of different types of 


active galaxies in this context is still missing. Here, we there¬ 
fore aim at a comparison of the molecular ISM in a sample of 
galaxies that includes starbursts, AGNs and ULIRGs. 

In this paper we present the first systematic unbiased spec¬ 
tral line survey of a sample of galaxies selected with respect 
to their nuclear activity. It is organised as follows; In Sect. [2] 
we introduce our galaxy sample. Section [3] describes the ob¬ 
servations and data reduction. Section [4] details the data anal¬ 
ysis, including a discussion about the molecular source sizes 
of the galaxies (Sect. 14.lt . the line identification and line pro¬ 
files (Sect. l4~2l >. spectral model fits, and rotational temperatures 
(Sect. l4~3l >. In Sect. 0 we summarise the detections (Sect. BTTl) . 
and present the column densities calculated under the LTE as¬ 
sumption (Sect. 15.2b . We indicate the potential key species for 
further extragalactic studies in Sect. [ 6 ] In Sect. [7] we detail our 
findings for starburst galaxies (Sect. mi). AGNs (Sect.im and 
ULIRGs (Sect. 17.3b . Finally, we summarise our main conclu¬ 
sions in Sect. [ 8 ] 


2. The galaxy sample 

We selected a sample of eight galaxies based on their luminosity 
and galactic nuclear activity, following these characteristics: (a) 
the galaxies are in the local universe (z < 0.05, D < 200 Mpc), 

(b) they show bright molecular emission (r™ (1_0) >50 mK) and 
cover a wide range of infrared luminosity (Lir = [10 9 -10 12 ] L e ), 
thus including low luminosity galaxies, LIRGs and ULIRGs; 

(c) the sample encompasses the two main types of galactic nu¬ 
clear activity: AGNs and starbursts; (d) it includes starbursts 
and ULIRGs in different stages of evolution; (e) it includes 
AGNs of both Seyfert (Sy) 1 and Sy 2 types. The eight galaxies 
are well-known extragalactic objects studied across almost the 
whole electromagnetic spectrum, and in particular at mm/sub- 
mm wavelengths. Moreover, these galaxies are often considered 
as archetypes in the local universe of the different nuclear activ¬ 
ities/stages of evolution mentioned above. Some of their main 
properties are shown in Table Q] 

Our sample includes the three starburst galaxies M83, 
NGC253, and M82. The last two were targets of line sur¬ 
veys at A = 2 and 1 . 3 mm done with the IRAM 3 0 m tele¬ 
scope bv IMartln et alJ (l2006bl) and lAladro et alJ (1201 lbl) . These 
two studies confirmed that the molecular ISM in the north¬ 
eastern molecular lobe of M82 shows typical features of an 
evolved starburst, where the gas is heavily pervaded by strong 
UV fields, while the nuclear region of NGC253 is more in¬ 
fluenced by lo w-velocity shocks. NGC 253 hosts nota ble large 
scale outflows dTurner & Hoi 119851: iBolatto et alJ 1201 3l) which 
indicate that its starburst is not in its initial state, though not 
as evolved as the M82 one. With the present A-3 mm obser¬ 


vations, these two galaxies are the best studied at millimetre 
wavelengths so far. We also observed M83 as an example of 
a young starburst (Mar tin e t al.l 2009c ). It h as a lower star for¬ 
mation rate (2.5 M-.~ vG 1 . 1 Walter et al.l 20081) than NGC 253 and 
M 82, and a star formation episode lasting (only) around 6 Myrs 
dHoughton & Thattel 120081) . Thus its chemical composition is 
expected to be somehow different from those of NGC 253 and 
M 82. With these three starburst galaxies, we aim to characterise 
better the changes of the molecular gas composition as the star- 
burst processes evolve. 

As galaxies containing AGNs, we include NGC 7469 (Sy 1), 
NGC 1068 (Sy 2), and M 51 (Sy2). The results of the NGC 1068 
survey were already published in lAladro et al.l ( 2013 ). and we re¬ 
fer the reader to that paper for a discussion of its isotopic line ra¬ 
tios and a chemical modelling of its central gas. Both NGC 7469 
and NGC 1068 have a starburst ring at ~lkpc from the cen¬ 
tral AGN, and have very similar observed properties (see e.g. 
Table 2 of IWilson et al.lll99ll) . On the other hand M51 shows 
sever al star-forming regions along its spiral arms (Bastian et al. 
12005' ). In the three cases, the angular resolution of the 30 m 
beam (~ 25") did not allow us to separate the emission coming 
from the central AGNs and those of the surrounding starbursts. 
Thus, our data might be showing not only the molecular com¬ 
position of their AGN but also some contribution from the star- 
burst regions. This might be particularly significative in the case 
of NGC 7469 due to its farther distance, although estimations of 
the emission coming from starburst regions at our frequencies 
are only available for NGC 1068 (see details in Sect. 17.21) . 

To complement our sample, we observed two of the clos¬ 
est ULIRGs, Arp 220 and Mrk231. Vigorous star forming re¬ 
gions and black holes coexist in the centres of ULIRGs. Still, 
the main powering source of their nuclei remains unclear in 
most cases. Arp 220 is in an early stage of the merging pro¬ 
cess, with a double nucleus that is not resolved by our obser¬ 
vations. An interferometric line survey done towards this galaxy 
bv lMartln et al.l (1201 ll) pointed to starburst processes as the more 
likely heating mechanism of its two centres (see also lGreve et al.l 
120091) . In contrast, Mrk231 is in a later stage of merging, where 
the two nuclei of the individual galaxies are already fused, and 
molecular-rich large-scale outflows and super-winds are clearly 
present (Aalto et al. 20121). AGN activity likely prevails i n this 
galaxy (e.g. lGallagher et all2002l:rvan der Werf et al.ll2010l). an d 
its nucleus seems to host a Sy 1 object (Boksenb erg et ail!l977h . 


3. Observations and data reduction 

The observations were carried out with the IRAM 30 m tele- 
scoptfl (Pico Veleta Observatory, Spain) between June 2009 
and March 2012. We observed the eight galaxies in the fre¬ 
quency range ~[86-116] GHz. For NGC 253 and M82 we also 
included lower frequencies from 80 GHz and 85 GHz, respec¬ 
tively. The reduced spectra are shown in Figs. Q]- 0 The half¬ 
power beam width (HPBW) corresponding to the survey fre¬ 
quencies ranged from_29" to 21". We used the band E0 of the 
EMIR receiver ( Carter et al.ll2012| ) and the WILMA autocorre¬ 
lator. This receiver-backend configuration allowed us to cover 
8 GHz simultaneously in the two orthogonal linear polarisations, 
and led to an original channel-width spacing of 5 - 7 km s^ 1 . The 
observed ~30GHz were covered by 5 tunings with a minimum 
overlap of 0.8 GHz. For some galaxies, an extra tuning has been 
used to increase the signal to noise in parts of the spectrum show- 

1 IRAM is supported by INSU/CNRS (France), MPG (Germany), 
and IGN (Spain). 
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ing faint lines. The observations were done by wobbling the sec¬ 
ondary mirror with a switching frequency of 0.5 Hz and a beam 
throw of 110" in azimuth (80" in the case of Arp 220). 

The weather was variable during the different observing 
runs, from excellent (pwv~lmm) to mediocre (pwv~6mm) 
weather conditions. We checked the pointing accuracy every 
hour towards several nearby bright continuum sources. The 
pointing corrections were always <5". The focus was also 
checked and corrected at the beginning of each run and during 
sunsets and sunrises. 

The data were first calibrated to the antenna temperature 
(Tn scale using the chopper-wheel method dPenzias & Burrusl 
119731) . The observed spectra were then converted to main beam 
temperatures (Tmb) using the relation Tmb = ( F e ff/B e g)T 
where F e ff is the forward efficiency of the telescope, whose val¬ 
ues were between 0.94 and 0.95, and B e ff is the main beam effi¬ 
ciency, ranging from 0.77 to 0.810. 

Each individual spectrum was analised separately. After 
eliminating bad channels or spectra containing strong ripples, all 
the spectra tuned at the same frequency were averaged and base¬ 
lines of orders between zero and two were subtracted to each 
chunk of 8 GHz. We note however that in the case of M 82 the 
baseline correction was not good enough to remove some small 
ripples in the range ~ [107 - 109] GHz. Higher order baselines 
did not improve the resulting spectra. Fortunately, comparing 
with the rest of the galaxies, no important lines are expected in 
that range. 

Due to the broad line-widths of these galaxies (>100 km s~ 1 ), 
the final spectra of M51, M82, M83, NGC 1068 and NGC253 
were smoothed to 11 - 14kms _1 , those of Arp 220 and 
NGC 7469 were smoothed to 22 - 27kms ', and the Mrk231 
data to ~60 km s 1 . The rms achieved for all galaxies was in av¬ 
erage 1 - 2mK (< lOmJybeam ') across the whole survey at 
the final velocity resolutions. The data were also corrected to 
first order for beam dilution effects as 


Th = + 0 2 b ) / eh Tmb , (i) 

where 7 b is the source averaged brightness temperature, 6 S is the 
molecular source size (see Table Q]and Sect. l4.H for details), and 
0b is the beam size. Gaussian distribution of the emission was 
assumed. We did not correct the data for the gain-elevation curve 
and for the gain losses due to a large wobbler throw because 
these corrections are only minor at 3 mm wavelength. 


4. Data analysis 

4. 1. Source sizes 

The source size of the molecular emission in each galaxy can 
be estimated in two ways. One is through high angular resolu¬ 
tion interferometric maps of the region of all the species, and the 
other is by comparison of the observations carried out with two 
(or more) telescopes having different beam sizes at the same fre¬ 
quency (so 0 S can be estimated from Eq. |T}. Unfortunately, the 
interferometric maps performed so far for our sample of galax¬ 
ies encompass a few abundant species only, such as CO, 13 CO, 
or HCN. Besides, most of the line transitions studied here were 
previously observed (if detected at all) with only one single-dish 
telescope. 

To simplify the situation, we adopt the same source 
size for all the molecules detected in a galaxy. We esti¬ 
mated their 0 S from interferometric maps of 12 CO(l - 0), 

2 http://www.iram.es/IRAMES/mainWiki/Iram30mEfficiencies 


HCN 

Ooo) 

1 - 0), HCO + (1 - 0), 13 CO(l - 0), and HCO(1 0 , - 

Downes & Solomon 1998 for Arp220 and Mrk231, 

Kohno et al. 

1996 for M51, Garcfa-Burillo et al][2002 for the 

NE lobe of M 82, Muraoka et al. 2009 for M 83, Heifer & Blitz 

1995, Schinnerer et al. 2000 for NGC 1068, Knudsen et al. 2007 

for NGC 253, and Davies et al. 2004 for NGC 7469). The values 


are listed in Table Q] 


The source size has a direct impact on the resulting bright¬ 
ness temperature (Eq. Q3, and consequently on the source aver¬ 
aged column density.The larger (smaller) the source size that is 
assumed, the lower (higher) the column densities that are de¬ 
rived. To work around this, we base our discussion on column 
density ratios between species. Our estimates of these ratios 
would be only sensitive to differences in the emission extent be¬ 
tween species, on which there is not enough information. 


4.2. Line identification , line profiles , and Gaussian fits 

We used the MADCUBAJJ s oftwarJl which i n cludes 
the CDMS and J PL catalogues dMuller it al.l 1200 ll 120051: 
iPickett et al.l Il998l) to identify the line s. The procedure 


similar to what is described in detail in Aladro et al 


and also shown in the NGC 1068 survey (lAladro et al 


ionS), 


20131) . 


We used CLASS0 to check whether there were lines com¬ 
ing from the image band of the spectra. The image band rejec¬ 
tion of the E0 receiver is >10 dB, which ensures that most of the 
lines from the image band, except for the strongest ones, are not 
entering the signal band. We found that only the spectra from 
NGC 253, M 83 and M 82 are affected by features coming from 
the image band. Few of them were blended with other transi¬ 
tions from the signal band. For those, whenever the separation 
was large enough, we fixed the observed line parameters in the 
signal band to estimate their contribution to the image band line. 
Details of the fitting parameters are shown in Appendix lA] 

Among the eight galaxies, only M83 and Mrk231 show 
simple Gaussian line profiles. The bright lines in Arp 220 (e.g. 
HCN, HCO + , CS) have a clear double-peak profile. It is not 
clear whether this might be due to the unresolved double nucleus 
or to absor ption effects, as reported in previous inte rferomet- 
ric studies (ISakamoto et al.l[20091 iTunnard et al.ll2015h . Fainter 
lines also show a deviation from Gaussian profiles, though not 
so clear as those which are optically thicker. Mrk 231 is also ex¬ 
pected to be affected by absorption, but its lines are compara¬ 
tively much fainter, and thus this effect, if present, is not so ev¬ 
ident. Regarding NGC 253, the observed circumnuclear region 
has five clumps of warm dense gas, but interferometric obser¬ 
vations are needed to resolve them (e.g. ISakamoto et ahll2011l) . 
With the 30 m telescope, the integrated spectrum of NGC 253 
appears as two blended Gaussians. On the other hand, the ‘two¬ 
horned’ line profiles of M51 and NGC7469 include not only 
the emission from the very centre, but possibly also that com- 
ing from the inner arm s and the surroundi ng molecular disks 
(Mat sushita et al] 1 19981: iDavies et alJ l2004l) . NGC 1068 has a 
similar physical structure as NGC 7469 (a circumnuclear disk, 
surrounded by a starburst ring at ~ 1 kpc from the centre). In ad¬ 
dition, the gas dynamics in the CND of NGC 1068 are quite com¬ 
plex, and some of the brightest lines, as seen by high-resolution 
observations, need sever al (two to four) Gaussian components 
to be tter fit the profiles (Krips et aljl201 It iGarcf a-Burill o et alJ 
120141) . Finally, while our M 82 observations are centred in the 
North-Eastern (NE) molecular lobe of the nucleus, some emis- 


3 http://cab.inta-csic.es/madcuba 

4 http://www.iram.fr/IRAMFR/GILDAS/ 
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Table 1: Main properties of the galaxies 


Galaxy 

a (J 2000) 
h:m:s 

(5 (J 2000) 

D" 

Mpc 

V b 

V LSR 

km s 1 

% 

"/pc 

SFR d 

M 0 yr 1 

^IR 

Lq 

Activity type 

M83 

13:37:00.93 

-29:51:56.40 

4.5 

513 

15/327 

2.5 

1.2 x 10 lu 

Starburst 

NGC 253 

00:47:33.12 

-25:17:18.60 

3.9 

250 

20/378 

3.6 

2.8 x 10 10 

Starburst 

M 82-NE 

09:55:51.90 

69:40:47.00 

3.5 

300 

12/204 

10 

5.9 x 10 10 

Starburst 

M51 

13:29:52.70 

47:11:43.00 

8.4 

470 

15/611 

2.5 

2.6 x 10‘° 

AGN 

NGC 1068 

02:42:40.90 

-00:00:46.00 

14.4 

1100 

4/279 

0.4 

1.9 x 10“ 

AGN + Starburst 

NGC 7469 

23:03:15.60 

08:52:26.00 

29.5 

4892 

6/858 

30 

3.9 x 10“ 

AGN + Starburst 

Arp 220 

15:34:57.08 

23:30:11.30 

70.0 

5350 

2/679 

240 

1.6 x 10 12 

Starburst 

Mrk 231 

12:56:14.20 

56:52:25.00 

170.0 

12173 

2 / 1640 

220 

3.2 x 10 12 

AGN 


For M 82 we observed the northeastern (NE) molecular lobe at an offset position (Aa-, AS) = (+13.0", +7.5") with respect to its dynamical 
centre. “ Distances taken from the NASA/IPAC Extragalactic Database; b Systemic Local Standard of Rest (LSR) radial velocities following 
the optical conven tion, taken from SIM BAD.S ource sizes and their co rresponding spatial scale s for the assumed dis tances. d St ar formation 
rate obtained from lWalteretalj j[200f ) for M 83: IStricklan d et all d2004|) fo r NGC 2 53 and M 82: ISchusteretaL[ j2007 ) for M 51: lEsquei et all 
d2014l) for NGC 1068: iGenzel et all d 19951) for NGC 7469: lAnantharamaiah et all 1200(f) for Am 220: iTavlor et all 119991 for Mrk231. We 


note that this parameter sho uld only 
luminosities taken from lSanders et al 


ive a rough idea of the activity, as it was calculated for different volumes in each galaxy. c Infrared 


give a r 


sion from the very centre enters the 30 m beam. This translates 
into a bump at low velocities of some of the lines that we can 
easily separate from the NE emission (see I Aladro et al.l [201 lbl 
for details). To homogenize the data, we fitted a single Gaussian 
to all lines detected in the survey. This simplification has little 
i mpact on the derived colu mn densities of the sp ecies, as shown 
bv lMartm et alJ (1201 ll) and lAladro et al.l (1201 lal) for Arp 220 and 
NGC 253 respectively. However, we made sure that the results 
obtained by fitting one or more Gaussians to the lines give dif¬ 
ferences in the overall column densities that are less than 10 % 
of their values. 


4.3. Spectral line model fit 

Molecular emission was modelled using the SLIM package 
within MADCUBAJJ (Martin et al. in prep.). The fit is per¬ 
formed in the parameter space of column density, rotational tem¬ 
perature, velocity and width of the line to the actual emission. 
Thus integrated and peak line intensity, as well as opacity, are 
derived from those parameters. The modelling assumes LTE, 
but not optically thin emission. However, the fixed values of the 
source sizes and rotational temperatures make the fitted lines op¬ 
tically thin and the line profiles Gaussian. A Gaussian fit has 
been performed to the CO(l - 0) lines, but as it is highly af¬ 
fected by opacity, which impact the column density estimate, 
we do not include CO in our results and discussion. For the rest 
of the species, we fitted the molecules one by one, including in 
each case all the transitions along the observed frequency range 
(in case there are more than one). As an exception, blended fea¬ 
tures were fitted simultaneously in order to take into account the 
contribution of each one to the total integrated areas. 

The rotational temperature ( T mt ) of a given species can only 
be constrained if there are observations of several rotational tran¬ 
sitions. Additionally, if the molecule is faint, the error associated 
to its calculated T rot could be, by far, larger than the value itself. 
As most of the species have only one transition in our survey, we 
fixed 7 , rot =10K. On average, as seen in our previous works, the 
rotational temperatures of the detected species are not expected 
to be much lower or higher than this value, though some excep- 


tions are observed, such as seen for NHiCN and CFLC’CH in 

NGC 253 and M82 (Martin et al. 

2006b; Aladro et a]J 20111), 

HC 3 N and C-C 3 H 2 in Arp 220 

Martin et al. 2011), and CN 


in Mrk231 and Arp 220 (iHenkel et al.l 12014b . The uncertainty 
added to the calculated column densities by assuming a fixed 
rotational temperature is not critical when comparing fractional 
abundances among galaxies; multi-line observations of many 
molecules in the nuclei of different galaxies and in the GC show 
that species have also very similar excitation conditions in differ- 
ent so urces (e.g. Aladr o et al.l 1201 lbl lArmiios-A bendan o et al.l 
12015D . Thus the associated rotational temperature of a given 
molecule is expected (or measured) to be similar for all the 
galaxies of our sample, being the variations of the T I0t propor¬ 
tional to the differences in the gas temperature among sources. 
For this reason, comparing fractional abundances instead of col¬ 
umn densities, makes the analysis less sensitive to these changes. 

Apart from the species mentioned before, there is the special 
case of the vibrationally excited lines of HC 3 N, detected (only) 
in Arp 220. They are expected to trace gas with tempera tures 
above 100K dCostagliola & Aahol2010tlMartm et al J2Q 1 ll) . For 
those vibrational transitions, we did calculate the temperature in 
a similar way as for the rotational temperatures, and obtained 
T vib = 190 +20 K. 

Warm kinetic temperatures (7f m ) estimated in the 
nuclear regions of these galaxies, such as 120 K for 
NGC253 (Bradford et al] [2003) or >150K for Mrk231 
(Ivan der Werf et al.l 20 l d) 7 come from higher frequency obser¬ 
vations and are representative of a small fraction of the gas 
only. Therefore, our assumption is still valid for the bulk of the 
molecular gas. Note, in any case, that rotational temperatures 
are lower than 7f m , as the H 2 densities are below the critical 
densities. Also, previous extragalactic surveys used similar 
assump tions, such as T rn i-1 2+6K for NGC253 dMartin et al.l 
H), 20 + 10 K for M82 (IXiadroetaD l 20 Tlbl) . and 10+5K 


assumpi 

12006hi) .__ - 

for NGC 1068 (lAladro et al.ll2013l) . 


5. Results 

5. 1. Summary of detections 

In ag reement with our previous surveys dAladro et al.l 1201 lbl 
Him we consider “detected” lines as those with a signal to 
noise ratio SNR > 3 and whose profiles can be clearly fitted by a 
Gaussian. Lines are considered “tentatively detected” when hav¬ 
ing 1 < SNR < 3, but can also be fitted with Gaussian profiles, 
whose velocities and line widths are in good agreement with the 
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ones obtained for detected lines (e.g. HC 5 N lines, see Fig. (3). 
“Non-detected” lines are those with SNR < 1 for which it is not 
possible to fit reliable Gaussian profiles, as the lines are domi¬ 
nated by the spectral noise. 

Thirty-seven species (including isotopologues) were iden¬ 
tified in NGC253, which make it the most prolific galaxy of 
our sample in terms of molecular complexity. We detected the 
13 C bearing isotopologues of CO, HCN, HCO + , HNC, CN, and 
CS, the ls O bearing isotopologues of CO and HCO + , the n O 
isotopologue of CO, and the 34 S bearing isotopologue of CS. 
Furthermore, 13 CN, S0 2 , NH 2 CN, HOCO + , HC 5 N, H 2 CS, and 
C 2 S were only detected in NGC 253. Tables[2]aiid Q]list the total 
column densities of the thirty-six molecules and isotopologues in 
the eight galaxies (excluding CO, see Sect. 14.21 ). Molecules de¬ 
tected (tentatively or not) for the first time in each of the galaxies 
are marked in boldface. 

Although tentative, this is the first time that HC 5 N is detected 
outside the Milky Way. All the seven brighter transitions based 
on our synthetic spectra of HC 5 N appear with measured inten¬ 
sities at 1 - 3 cr level (see Table [A2l and Fig. [3- Deeper ob¬ 
servations would be needed to confirm this. We note, however, 
some inconsistencies in the identification of this species between 
iMeier et al . (1201 5) and this work. The features that we identify 
as the HC 5 N transitions (J-J')-( 33-31), (37-36), (38-37) at 
87.8, 98.5, and 10L2GHz respectively, are tentatively identified 
by Meier et al . ( 2015h as NH 2 CHO, CH 3 CH 2 CN, and CH 3 SH, 
also in the centre of NGC 253. We do not detect emission of 
those molecules in our survey, although we cannot rule out that 
this might be due to differences in sensitivity between the IRAM 
30 m telescope and ALMA. However, we cover a broader range 
of frequencies where other HC 5 N lines are also consistently de¬ 
tected, which makes us favour our identification. HC 5 N was pre¬ 
viously observed towards Galactic dense cores of dark clouds, 
and seems to trace the initial conditions of star formation (e.g. 
ISuzuki et al.lll992l : lRathborne et al.ll2008l) . 

In addition, we found that c-C 3 H 2 , NS and, especially, 
CH 3 CHO were detected in NGC 253 and M 82, but not in the 
rest of the galaxies of our sample. Outside the Milky Way, 
CH 3 CHO was also tentatively detected in NGC 253 (iMeier et al. 
2015 ) and towards a galaxy at redshift z=0.89 dMuller et al. 
2011). Although not much is known about this molecule, 
its origin seems to be linked to grain mantle destruction by 
shocks, and is found in massive star-formin g regions (e.g. 
IChengalur & Kanekarll2003t iBennett et al.l[2005l) . 

As commented in Sect. 14.31 we detected vibrationally excited 
lines of HC 3 N only in Arp 220, which are listed in Table IA.7I 
They trace the hottest and densest molecular gas in the inner 
regions of hot cores in star-form ing regions (Ide Vicente et al.1 
2000; Mar tin-Pintado et akll2005l) . 

5.2. Column densities and fractional abundances 

For those species not detected in a galaxy, we estimated the 3<x 
upper limit to the column density by assuming a limit to the in¬ 
tegrated area (I) given by 

I < 3 x nns x V FWHM X D, (2) 

where rms is the lcr noise level of the spectrum that contains 
the line, Z) v is the spectral resolution in velocity units indicated 
in Sect.[3for each galaxy, and FWHM is the average full width 
at half maximum of the Gaussian fits measured for other species 
in the same galaxy. As FWHM we generally used 95 kms -1 for 
M 83, 200km s - 1 for NGC 253, 100km s ' 1 for M 82, 130km s ' 1 


for M51, 240kins' 1 for NGC 1068, 250kms” 1 for NGC7469, 
390 kms ^ 1 for Arp 220, and210kms~' for Mrk 231. Note, how¬ 
ever, that these values may differ for some species. 

We calculated fractional abundances of all the species (or 
their upper limits in case they were not detected) relative to the 
12 C 18 0 (hereafter C 18 0) column density, which is significantly 
less affected by opacity effects than the main oxygen isotopo¬ 
logue I 2 C ie O. We note, however, that C ls O seems to be en¬ 
hanced in the two ULIRGs of our sample. In Sect. 17.31 we dis¬ 
cuss how this might affect our comparison of abundances among 
galaxies. 


6. Key molecular species 

Figure [3 shows the comparison of thirty-five molecular abun¬ 
dances among the eight galaxies in the sample. Table |4] lists the 
most relevant column density ratios used for the discussion in 
this section and the following. Based on these results, we indi¬ 
cate here the species presenting the highest abundance contrasts 
among our galaxy sample. These molecules are potentially good 
candidates for further high-resolution follow up observations. 
We note that we only discuss here species where the abundance 
difference is at least a factor of three. 


- CH 3 CCH is detected in M 83, NGC 253, M 82, and Arp 220, 
but not in the AGNs of our sample. M82 has the max¬ 
imum CH 3 CCH abundance and, in fact, no other galaxy 
has been found to have such a high value (a t l east re lative 
to H 2 , C 34 S and C ls O, lAladro et al l [201 1 alibi [2011) . The 
other galaxies where this species was detected have sim¬ 
ilar values, five to seven times lower than that of M82. 
Interestingly, CH 3 CCH was not detected in NGC 1068, 
NGC 7469, M51, and Mrk 231, which supports the idea that 
this molecule is highly contrasted between starbursts and 


gniy con 

AGNs (IAladro et all 12013 ). The upper limits to the abun¬ 
dance of M51 and NGC 1068 are significant, and show that 
the difference in abundance between M 82 and NGC 1068 is 
about one order of magnitude. 

In the Galaxy, CH 3 CCH was found to be very abundant 
in regions of prolific star formation, where its column 
density is substantially higher (MrmrrH ~ 10 14 ~ 15 cm 2 , 

in 


density is substantially mgner (1Vch,cch ~ 1 U cm 
IChurchwell & Hollisi 11983c iMiettinen etal. 2006j), than 


dark clouds, (Nch,cc h ~ 10 13 cm 2 , Irvine etaD 11981b 
iMarkwick et akll2005l) . Our values for M 82, M 83, NGC 253 
and Arp 220 are of the same order as the ones cited above 
for high-mass star-forming cores. However, given the differ¬ 
ent scales in extragalactic and galactic sources, the compar¬ 
ison of the column densities is not straightforward. In ad¬ 
dition, it was recently found that CH 3 CCH is enhanced in 
PDRs in the Horsehead nebula as compared to dense cores 
( iGuzman et al.l 120141) . High temperatures could help to un¬ 
derstand high abundances of CH 3 CCH in these objects, as 
it may form after the desorption of its chemical precur¬ 
sors from grain mantles. In fact, rotational temperatures de¬ 
rived for Orion-KL, SgrB2, DR 21 and many other simi¬ 
lar sources are in the range [20-50] K (IChurchwell & Hollisi 

1 19831 IMiettinen et al.ll2006l). _ 

As previously claimed by lAladro et all (2013), and in agree¬ 
ment with our results here, the abundance of CH 3 CCH seems 
to be enhanced in starburst galaxies. It is a known tracer of 
relatively dense gas (the critical densities for the two tran¬ 
sitions observed here are 10 4 - 10 5 cm -3 ). Its emission is 
more widespread than that of other similar molecules, such 
as CH 3 CN, and can also be found in regions outside the cores 
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Aladro et al. (2015): A =3 mm line survey of nearby active galaxies 
Table 2: Source averaged column densities [cm" 2 ] 


Molecule 

M83 

NGC 253 

M 82 

M51 

rj CO 

(5.7 ± 0.1) X 

10 lb 

(1.88 ± 0.03) x 

10 u 

(1.32 ± 0.02) x 

10" 

(2.1 ± 0.1) x 10 lb 

C 18 0 

(1.10 ± 0.03) x 

10 16 

(5.1 ± 0.2) x 

10 16 

(2.61 ± 0.08) x 

10 16 

(5.2 ± 0.4) x 10 15 

C 17 0 

(1.9 ± 0.2) x 

10 15 

(6.7 ± 0.3) x 

10 15 

(1.6 + 0.3) X 

10 15 

(9.0 ± 1.9) x 10 14 

HCN 

(8.6 ± 0.1) x 

10 13 

(3.85 ± 0.05) x 

10 14 

(2.90 ± 0.06) x 

10 14 

(4.5 ± 0.1) x 10 13 

h 13 cn 

( 2.3 ± 0 . 5 ) x 

10 12 

(2.8 ± 0.1) x 

10 13 

(5.0 ± 0.6) x 

10 12 

(1.8 ± 0.4) x 10 12 

HNC 

(3.13 ± 0.09) x 

10 13 

(1.81 ± 0.02) x 

10 14 

(1.16 ±0.03) x 

10 14 

(1.35 ± 0.06) x 10 13 

HN 13 C 

( 1.2 ± 0 . 6 ) x 

10 12 

(1.5 ± 0.1) x 

10 13 

< 2.1 x 

10 12 

< 1.4 x 10 12 

HCO + 

(4.6 ± 0.1) x 

10 13 

(1.94 ± 0.03) x 

10 14 

(2.55 ± 0.06) x 

10 14 

(1.31 ±0.04) x 10 13 

H 13 CO + 

( 1.1 ± 0 . 4 ) x 

10 12 

(1.3 ± 0.2) x 

10 13 

(6.9 ± 0.3) x 

10 12 

< 6.6 x 10“ 

hc 18 o + 


- 

(4.3 ± 0.2) x 

10 12 

(4.7 ± 2.5) x 

10 12 

- 

HOC + 

( 1.3 ± 0 . 2 ) x 

10 12 

(6.4 ± 0.4) x 

10 12 

(3.6 ± 0.2) x 

10 12 

< 5.8 x 10“ 

HCO 

< 5.0 x 

10 12 

(5.9 ± 0.2) x 

10 13 

(5.1 ± 0.6) x 

10 13 

< 7.0 x 10 12 

c 2 h 

( 5.0 ± 0 . 1 ) x 

10 14 

(2.82 ± 0.04) X 

10 15 

(3.69 ± 0.05) x 

10' 5 

(1.6 ± 0.1) x 10 14 

CN 

( 3.24 ± 0 . 04 ) x 

10 14 

(1.74 ± 0.01) X 

10 15 

(1.02 ± 0.02) x 

10 15 

(1.31 ±0.02) x 10 14 

13 CN 

< 7.0 x 

10 12 

(6.2 ± 0.5) x 

10 13 

< 1.2 x 

10 13 

< 1.5 x 10 13 

ch 3 oh 

( 1.9 ± 0 . 2 ) x 

10 14 

(1.3 ± 0.1) x 

10 15 

(2.07 ± 0.08) x 

10 14 

(3.9 ± 1.5) x 10 13 

cs 

(5.7 ± 0.1) x 

10 13 

(4.07 ± 0.07) X 

10 14 

(2.50 ± 0.05) x 

10 14 

(1.8 ± 0.1) x 10 13 

c 34 s 

(7.7 ± 1.0) x 

10 12 

(5.4 ± 0.2) x 

10 13 

(1.50 ± 0.08) x 

10 13 

(3.1 ± 0.8) x 10 12 

13 cs 

< 2.9 x 

10 12 

(1.4 ± 0.1) x 

10 13 

< 4.6 x 

10 12 

< 3.2 x 10 12 

hc 3 n 

( 1.9 ± 1 . 2 ) x 

10 13 

(1.04 ± 0.05) x 

10 14 

(4.6 ± 0.4) x 

10 13 

< 4.4 x 10 12 

so 

( 2.0 ± 0 . 3 ) x 

10 13 

(1.7 ± 0.3) X 

10 14 

(7.7 ± 0.8) x 

10 13 

(1.1 ± 0.3) x 10 13 

so 2 

< 1.6 x 

10 13 

(3.6 ± 0.6) x 

10 13 

< 2.0 x 

10 13 

< 1.3 x 10 13 

NS 

< 1.2 x 

10 13 

(3.8 ± 1.9) x 

10 13 

(3.0 ± 0.6) x 

10 13 

< 6.4 x 10 12 

HNCO 

(4.7 ± 2.5) x 

10 13 

(2.2 ± 0.7) x 

10 14 

( 2.7 ± 5 . 0 ) x 

10 13 

< 6.7 x 10 12 

n 2 h + 

( 6.5 ± 0 . 3 ) x 

10 12 

(3.97 ± 0.07) x 

10 13 

(1.31 ± 0.03) x 

10 13 

(4.1 ± 0.3) x 10 12 

SiO 

( 2.4 ± 0 . 6 ) x 

10 12 

(3.1 ± 0.3) X 

10 13 

(4.4 ± 0.6) x 

10 12 

< 1.2 x 10 12 

CH 3 CN 

( 3.4 ± 0 . 9 ) x 

10 12 

(3.3 ± 2.6) x 

10 13 

(4.8 ± 15.8) x 

10 12 

< 1.8 x 10 12 

ch 3 cch 

( 9.4 ± 5 . 7 ) x 

10 13 

(5.4 ± 1.0) x 

10 14 

(1.15 ± 0.04) x 

10' 5 

< 5.3 x 10 13 

c-C 3 H 2 

< 3.7 x 

10 15 

(1.05 ± 0.06) x 

10 14 

(1.31 ±0.06) x 

10 14 

< 1.7 x 10 15 

nh 2 cn 

< 3.2 x 

10 12 

(1.3 ± 0.3) x 

10 13 

< 3.9 x 

10 12 

< 4.0 x 10 12 

HOCO + 

< 6.4 x 

10 12 

(1.2 ± 0.7) x 

10 13 

< 6.1 x 

10 12 

< 5.0 x 10 12 

OCS 

< 7.3 x 

10 13 

(4.9 ± 0.3) x 

10 14 

< 8.4 x 

10 13 

< 7.4 x 10 13 

C 2 S 

< 7.6 x 

10 12 

(2.5 ± 0.8) x 

10 13 

< 1.2 x 

10 13 

< 1.9 x 10 13 

H 2 CS 

< 1.5 x 

10 13 

(3.8 ± 0.9) x 

10 13 

< 1.4 x 

10 13 

< 1.1 x 10 13 

hc 5 n 

< 3.5 x 

10 14 

( 7.1 ± 13 . 4 ) x 

10 13 

< 3.8 x 

10 14 

< 4.0 x 10 14 

ch 3 cho 

< 1.1 X 

10 13 

( 7.4 ± 64 . 7 ) x 

10 13 

( 3.4 ± 0 . 4 ) x 

10 13 

< 1.5 x 10 13 


Molecules in boldface are new detections for a given galaxy. Column densities have been calculated for molecules clearly detected, as well as 
for those tentatively detected. Upper limits are shown for species not detected. 


of the molecular clouds (Ala dro et al.ll201 lah . The response 
of CH iCCH to the presence of UV fields and shocks has not 
been studied in detail in extragalactic sources. Consequently, 
the chemistry of CH 3 CCH is still uncertain, but it might 
well be that the paths for its production are more effective in 
the presence of UV radiation fields created by massive star 
formation (see I Aladro et al.lf201 lbi and references therein). 
Supporting this idea, recent interferometric maps of several 
species in the central half kiloparsec of NGC 253 show that 
the emission of CH 3 CCH h as a s imilar morphology to other 
PDR tracers, such as C 2 H dMeier et al.ll2015l) . Thus, it seems 
that, apart from being a dense gas tracer, CH 3 CCH can also 
be considered a PDR tracer. 

HC 3 N and H 13 CN show clearly enhanced abundances in 
Arp 220 and Mrk231 with respect to the rest of galaxies 
in our sample. This could be related to larger amounts of 
dense gas and warm dust present in L’LlRGs (as compared 
to less obscured galaxies. iLindberg et al.ll201 ll) . where the 
dust would prevent their photo-dissociation. On average, the 
differences between the ULIRGs and the other galaxies is a 
factor of three for HC 3 N and seven for H 13 CN, although the 


maximum differences between individual galaxies are of fac¬ 
tors four and twenty-one for HC 3 N and H 13 CN respectively. 
CH 3 CN and SiO have maximum abundances in Arp 220, 
which are a factor of ~7 higher than those of M 82, which 
shows the minimum in both cases. Both species are found 
to be tracers of dense gas and shocks, although enhanced 
SiO abun dances may also be rela t ed to X-ray chemis try in 
AGNs (iGarcfa-Burillo et aD 120101 : lAladro et al.ll2013h . The 
species are photo-dissociated in the external layers of molec¬ 
ular clouds if strong UV fields are present. This scenario is 
in agreement with both molecules having high abundances in 
Arp 220, which probably has more dense gas and dust, and a 
minimum in M 82, as it is heavily pervaded by UV fields. 
C 2 H. Previous observational studies found high abun- 
danc es of this species in PDRs of starburst galaxies 
(e.g. Meier & Turnedl2005t iFuente et~al]l2005b lAladro et~akl 
1201 lbl) . Together with M 82, the highest abundances of C 2 H 
in our sample are also those of NGC 7469 and NGC 1068, 
while M 51 presents the lowest value, followed by Arp 220 
and M83. There is a factor five of difference between the 
abundances of M82 and M51. In the case of NGC 1068 
and NGC 7469, significant emission is likely to be arising 

























Aladro et al. (2015): A =3 mm line survey of nearby active galaxies 
Table 3: Source averaged column densities [cm" 2 ] 


Molecule 

NGC 1068 

NGC7469 

Arp 220 

Mrk 231 

n CO 

(4.6 ± 0.1) x 

10" 

(3.6 ± 0.2) x 

10 lb 

(6.6 ± 0.4) x TO 17 

(5.6 ± 0.9) x 10 16 

c I8 o 

(1.29 ± 0.04) x 

10 17 

(6.4 ± 1.3) x 

10 15 

(6.2 ± 0.4) x 10 17 

(4.9 ± 0.9) x 10 16 

c 17 o 

< 1.0 x 

10 16 

< 3.9 x 

10 15 

< 2.9 x 10 16 

< 2.2 x 10 16 

HCN 

(2.45 ± 0.04) x 

10 15 

(1.04 ± 0.08) x 

10 14 

(4.6 ± 0.1) x 10 15 

(8.5 ± 0.7) x 10 14 

h 13 cn 

(8.1 ±0.8) x 

10 13 

< 1.4 x 

10 13 

( 7.1 ± 0 . 5 ) x 10 14 

(2.0 ± 0.5) x 10 14 

HNC 

(7.0 + 0.1) x 

10 14 

(4.2 ± 0.5) x 

10 13 

(3.2 ± 0.2) x 10 15 

(2.5 ± 0.7) x 10 14 

hn 13 c 

(2.4 ± 0.6) x 

10 13 

< 1.6 x 

10 13 

< 1.3 x 10 14 

< 1.8 x 10 14 

HCO + 

(9.1 ±0.2) x 

10 14 

(7.4 ± 0.5) x 

10 13 

(1.12 ± 0.06) x 10 15 

(3.3 ± 0.3) x 10 14 

H 13 CO + 

< 9.0 x 

10 12 

< 8.0 x 

10 12 

(9.8 ± 3.2) x 10 13 

< 9.1 x 10 13 

HC ls O + 


- 


- 

- 

- 

HOC + 

(1.6 ± 0.4) x 

10 13 

< 6.8 x 

10 12 

< 5.6 x 10 13 

< 6.2 x 10 13 

HCO 

(2.6 ± 0.8) x 

10 14 

< 7.1 x 

10 13 

< 5.7 x 10 14 

< 8.0 x 10 14 

C 2 H 

(1.10 ±0.02) x 

10 16 

(1.0 ± 0.1) x 

10 15 

( 2.2 ± 0 . 1 ) x 10 16 

(3.4 ± 0.7) x 10 15 

CN 

(6.85 ± 0.06) x 

10 15 

(4.6 ± 0.2) x 

10 14 

(6.1 ±0.2)x 10 15 

(1.8 ± 0.1) x 10 15 

13 CN 

< 1.5 x 

10 14 

< 6.6 x 

10 13 

< 5.8 x 10 14 

< 1.4 x 10 14 

CH 3 OH 

(2.4 ± 0.2) x 

10 15 

< 1.8 x 

10 14 

(5.6 ± 0.6) x 10 15 

< 2.6 x 10 15 

CS 

(7.6 ± 0.2) x 

10 14 

(8.1 ± 1.0) x 

10 13 

(3.36 ± 0.08) x 10 15 

(3.5 ± 0.8) x 10 14 

c 34 s 

(1.3 ± 0.2) x 

10 14 

< 3.0 x 

10 13 

(9.6 ± 1.3) x 10 14 

< 2.3 x 10 14 

13 cs 

< 3.6 x 

10 13 

< 3.0 x 

10 13 

< 2.0 x 10 14 

< 2.0 x 10 14 

hc 3 n 

(1.50 ±0.07) x 

10 14 

< 5.0 x 

10 13 

(2.7 ± 0.2) x 10 15 

(2.1 ±0.5)x 10 14 

so 

(2.2 ± 0.4) x 

10' 4 

< 6.9 x 

10 13 

< 4.4 x 10 14 

< 5.0 x 10 14 

so 2 

< 7.4 x 

10 15 

< 1.1 X 

10 14 

< 9.2 x 10 14 

< 1.1 x 10 15 

NS 

< 9.6 x 

10 13 

< 8.0 x 

10 13 

< 3.1 x 10 14 

< 5.5 x 10 14 

HNCO 

(3.2 ± 0.4) x 

10 14 

< 6.0 x 

10 13 

(1.8 ± 0.9) x 10 15 

< 3.0 x 10 14 

n 2 h + 

(1.2 ± 0.5) x 

10 14 

< 8.2 x 

10 12 

(5.9 ± 0.3) x 10 14 

< 5.4 x 10 13 

sib 

(6.8 ± 0.8) x 

10 13 

< 1.5 x 

10 13 

(6.3 ± 0.7) x 10 14 

< 1.7 x 10 14 

CH 3 CN 

(1.0 ± 2.9) x 

10 14 

< 1.7 x 

10 13 

(9.2 ± 1.7) x 10 14 

< 9.1 x 10 13 

CH 3 CCH 

< 5.7 x 

10 14 

< 3.7 x 

10 14 

(3.7 ± 0.8) x 10 15 

< 4.0 x 10 15 

c-C 3 H 2 

< 2.7 x 

10 16 

< 3.9 x 

10 15 

< 4.4 x 10 16 

< 5.3 x 10 14 

NH 2 CN 

< 2.6 x 

10 13 

< 2.2 x 

10 13 

< 1.5 x 10 14 

< 1.8 x 10 14 

HOCO+ 

< 7.7 x 

10 13 

< 4.4 x 

10 13 

< 3.7 x 10 14 

< 3.5 x 10 14 

OCS 

< 8.9 x 

10 14 

< 8.2 x 

10 14 

< 6.1 x 10 15 

< 9.2 x 10 15 

C 2 S 

< 1.2 x 

10 14 

< 9.3 x 

10 13 

< 6.5 x 10 14 

< 7.1 x 10 14 

H 2 CS 

< 1.6 x 

10 14 

< 9.0 x 

10 13 

< 6.6 x 10 14 

< 7.6 x 10 14 

HCjN 

< 5.4 x 

10 15 

< 4.9 x 

10 15 

< 4.0 x 10 IS 

< 3.9 x 10 16 

CH3CHO 

< 4.6 x 

10 14 

< 1.2 x 

10 14 

< 7.2 x 10 14 

< 1.3 x 10 15 


Molecules in boldface are new detections in a given galaxy. Column densities have been calculated for molecules clearly detected, as well as 
for those tentatively detected. Upper limits are shown for species not detected. 


from the circumnuclear star forming rings, as in NGC 1097 
(IMartrn et al.ll2015l) . 

CN and HCN abundances stand out in NGC 7469, 
NGC 1068, Mrk231, and M 82, being the difference between 
these galaxies and the rest of the sample larger in HCN. In 
agreement w ith previous observations (e.g. Aalto et al. 2002; 
Aaltoll2008l) and model predictions dLepp & Dalgarno 1996; 


Sternberg & Dalgarnol 19951 iMeiierink & Snaaiis 1200.51) . our 

data indicate that CN and HCN could be effectively pro¬ 
duced in the X-ray and UV-field dominated regions expected 
in AGN and starburst galaxies. 

HNCO is tentatively detected for the first time in M 82. This 
species is claimed to be a good indicator of the evolution of 
nuclear starbursts in galaxies, with abundances decreasing 
toward the last evolutionary stages of the starbursts (in those 
late phases HNCO should be efficiently destroyed by means 
of reactions involving the abundant UV photons and/or ions). 
HNCO remained undetected in previou s studies of M82 
dMartm et al.ll2009al : I Aladro et al.ll201 lbl) . Our new tentative 
HNCO line detections in M82 have intensities of ~3mK, 
which are higher than the u pper lim it of <1.3 mK previously 
estimated by Martin et~al] (12009al) . Our HNCO abundances 


in M 82 are four times lower than in NGC 253 and M 83. This 
contrast is smaller than that found by IMartrn etal . ( l2009a ) 
for HNCO/C 34 S, which shows differences of nearly two or¬ 
ders of magnitude, but goes in the same direction. Therefore, 
our observations agree with HNCO being a good tracer of the 
evolutionary stage of starburst galaxies. 

- C 32 S/C 34 S. Both 32 S and 34 S are believed to be a prod¬ 
uct of hydrostatic or explosive oxygen-burning in massive 
stars and Type la supernovae (see [Chin et al.ll 19961 and ref¬ 
erences therein). However, not much is known apart from 
some Galactic observations which indicate that_ 32 S/ 34 S in¬ 
creases with the galactocentric radius dChin et al..' 1996 ). In 
external galaxies, more data is needed to know how this ratio 
varies from source to source. 

We found that the C 32 S (2 - 1)/C 34 S(2 - 1) line ratio is 
five times higher in M 82 (with the maximum value of of 
21.5 in our sample) than in Arp 220 (minimum value of 4.5, 
similar to that derived with SO bv lMartm et al.ll2011 ). This 
points to a higher CS opacity in Arp 220 where there is much 
more molecular gas. As a reference value, C 32 S/C 34 S in the 
Galactic centre (GC), the local ISM, and the Solar system 
is ~22, dFrerking et al.ll 1 980l) . M82 has a similar value to 
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the GC, which is ~3 times higher than that of NGC253 
(7.6 ). More soph i sticate d estimates by iMartfn et alJ (1201 Oh 
and lHenkel et alJ (l2014l> give values close to the local inter¬ 
stellar value for NGC 253 (>16). 

NGC 1068, M51, NGC253, and M83 have C 32 S/C 34 S J = 
2 —>1 line ratios of 6 - 8 , well below the GC and M 82, and 
closer to the value in Arp 220. Though these ratios point out 
towards the degree of processing of the ISM in the galactic 
nuclei, the values should be taken with care, as opacity on the 
main i sotopologue might play an important r ole in these re¬ 
gions (IMartm et alJl20ld : iHenkel et alJl2014l and references 
therein). Opacity becomes less of an issue for LIRGs, based 
on the 13 C0/C 18 0 sample measured by Costagliola et al. 

d2onh . 


7. Chemical caracterisation of our sample 

7. 1. Starburst galaxies 

The molecular gas in the nuclear region of M 83 is rich in some 
shock tracers and species formed in dust grains, such as HNCO 
and CH 3 OH (Fig.[ 6 ]and Table[4}, that are easily released through 
grain mantle evaporation by shock waves. Although the abun¬ 
dances of these two species are high relative to C ls O (10 ~ 3 
and 10 -2 for HNCO and CH 3 OH respectively), other shock and 
dense gas tracers are just marginally detected (e.g. SiO, CH 3 CN, 
and SO). No Hu recombination lines were detected in M 83, un¬ 
like in the more active starbursts (see SFR in Table Q} M 82 and 
NGC 253, in spite of all three galaxies having similar beam fill¬ 
ing factors. This suggests that big and medium-size HII regions 
and OB star clusters are relatively scarce in M83. Accordingly, 
the PDR-tracers C-C 3 H 1 and HCO were not detected. 

As for NGC 253, groups of OB stars create photo-dissociated 
regions in the outskirts of the molecular clouds. This is 
confirmed by previous detections of PDR tracers, such as 
HOC + , CO + , HCO, and c-C 3 H 2 ( Martin et al] 2009a ). which 
are easily formed by ion-molecule reactions in g as phase (e.g. 
iHollis & Church weill 1983llThaddeus et all 1985ft . However, the 
gas in NGC 253 is mainly characterised by the presence of many 
shock tracers and species formed in the icy mantles of dust 
grains. HNCO, CH 3 OH, and SO are a factor of two more abun¬ 
dant than in M 83. Also, some other faint shock tracers have been 
only detected in this galaxy, such as HOCO + and NS,_as well as 
the sulphu r bearing species S0 2 . H 2 CS, and C 2 S (Martin et al. 
120031 [20051 ) which got confirmed here (see Fig.| 6 ]l. 

The gas in the NE lobe of M 82 is strongly affected by UV 
fields which shape its chemical composition (e.g. 1 Aladro et al.l 
l2011bl). Molecular clouds are highly ionized and fragmented 
dFuente et al.ll2008l) and, according to our results, the PDR trac¬ 
ers c-C 3 H 2 , HCO, and HOC + have slightly higher abundances 
than in NGC 253 (but just a factor of two or less. See Table |4}. 
The species that shows the largest difference between M 82 and 
the other two starburst galaxies is CH 3 CCH, being in the for¬ 
mer five times more abundant (see Sect. [ 6 ] for more details). 
Still, some shock tracers, such as CH 3 OH, SiO, and HNCO, 
are detectable in the cores of the molecular clouds that are well 
shielded from the UV radiation (IMartm et alj|2006al) . However, 
their abundances are four times lower than in NGC 253, indi¬ 
cating that, although present, shocks do not represent the main 
heating mechanism of the molecular gas as compared to UV ra¬ 
diation. 


7.1.1. A possible chemical scenario of starburst evolution 

As commented in Sect. [2] M 83, NGC 253, and M 82 can be 
taken as examples of starbursts in a early, intermediate, and 
late stages of evolution. If so, it may be possible to link 
their evolutionary stages to the chemical composition of their 
molecular gas. This evolutionary s cenario has been progres¬ 
sively fine-tuned by several authors (IGarcfa-Burillo et al.ll2002l : 
Mart in et al .1 l2006ra lAladroetal.ll2011bl) . An early starburst is 
expected to be chemically characterised by the absence of most 
PDR tracers (HCO, C-C 3 H 2 ) and Her recombination lines, as well 
as the presence of only the most common shock-tracers (HNCO, 
CH 3 OH, SiO). Physically, first shocks could be caused by close 
encounters or mergers between galaxies, or collisions of giant 
molecular cloud complexes (for example, following the bar po¬ 
tential toward the centre of the galaxy. iBinnev et aU 199 ll) . 

In an intermediate starburst stage, with similar characteris¬ 
tics to NGC 253, the stellar population would contain more re¬ 
cently formed young massive stars and supernovae events that 
ionise the outer layers of the molecular clouds, as seen by the 
detection of the PDR tracers C-C 3 H 1 , HCO and HOC + . During 
this phase, PDR tracer abundance enhancements are still smaller 
than those caused by shocks, which would reach their maxi¬ 
mum influence over the ISM at this point. Molecules such as 
CH 3 OH, HNCO, and SO have higher abundances than in a pre¬ 
vious (and later) phase. Other shock tracers would be also de¬ 
tectable, specifically HOCO + and some sulphur bearin g sp ecies, 
such as H 2 CS, OCS, SO, S0 2 and C 2 S (see IMartm et al.ll2003l 
120051 and Fig.| 6 ]>, whose abundances are enhanced during the last 
evolutionary stages of hot cores. 

Finally, an advanced starburst such as that found in M82, 
would be characterised by super-winds created by massive stars 
and supernovae, as well as a larger population of high-mass X- 
ray binaries, which provide a modest X-ray emission, needed to 
under stand the abundance of some molecular species, such as 
CO + (ISpaans & Meiierinkll2007l: IChiang & Kondl201 ll) . Stellar 
winds enrich the ISM with both PDR tracers (HCO, HOC + , c- 
C 3 H 1 , and CH 3 CCH) and shock-tracers (CH 3 OH, SiO, OCS, 
SO), the first ones being more abundant, since UV fields are 
more intense than in previous stages. Less abundant, or more 
complex, shock tracers and sulphur species would be easily 
photo-dissociated by the radiation fields, and thus would become 
less abundant than in previous starburst phases or not detectable 
at all. 

7.2. Seyfert galaxies 

Among the galaxies in our sample, M51, NGC 1068, and 
NGC 7469 host an A GN (the first two are classified as Sy2, 
and the latter as Sv 1. lOsterbrock & Martel 199 3). A detailed 
study devoted to our results in NGC 1068 has been presented 
by lAladro et al.l ( 2013 ). The CNDs of these galaxies are sur¬ 
rounded by arms of active star forming regions at ~ 1-2 kpc 
distance from their centres. Given the angular resolution of the 
IRAM 30 m telescope at 3 mm wavelengths (22" -29"), our data 
embrace the CNDs as well as the surrounding regions, so some 
contribution from the star forming rings is expected in our data. 
Other observations of NGC 1068 indicate that some molecules 
show significant emission from the starburst ring (< 2 5% for 
SiO (2 - 1) and < 30% for HCO + (l - 0). iUsero et al.ll200l . 
while high transitions of the dense gas tracers HCN, HCO' and 
CS ar e hardly detected outside the CND ( Garcfa-Burillo ~etaP 
120141) . Recently, Takajima et al (in prep.) carried out a similar 
molecular line survey of NGC 1068 with the Nobeyama 45 m 


























































Aladro et al. (2015): A= 3 mm line survey of nearby active galaxies 
Table 4: Relevant column density ratios used for the discussion. 


Molecular ratio 

M 83 

NGC 253 

M 82 

M51 

NGC 1068 

NGC7469 

Arp 220 

Mrk231 

HNCO / C l8 0 

4.3 x ltT 8 

4.3 x 1(T 3 

1.0 x 10 -J 

< 1.3 x 10 

2.5 x 10“ J 

< 9.4 x 10~ 3 

2.9 x 10 -J 

< 6.1 x 10' 3 

ch 3 oh/c I8 o 

1.7 x 10~ 2 

2.5 x 1(T 2 

7.9 x 10“ 3 

7.5 x 10~ 3 

1.9 x 10~ 2 

< 2.8 x 10~ 2 

9.0 x 10“ 3 

< 5.3 x 10' 2 

SiO / C ls O 

2.2 x 1(T 4 

6.1 x 1(T 4 

1.7 x 10“ 4 

< 5.0 x 10~ 5 

5.3 x 10- 3 

< 2.3 x 10“ 3 

1.1 x 10~ 3 

< 3.5 x 10~ 3 

ch 3 cn/c 18 o 

3.1 x 1(T 4 

6.5 x 1(T 4 

1.8 x 10~ 4 

< 1.8 x 10~ 3 

7.8 x 10- 4 

< 2.6 x 10- 3 

1.5 x 10~ 3 

< 1.8 x 10~ 3 

ch 3 cch/c 18 o 

8.5 x 10- 3 

1.0 x 1(T 2 

4.4 x 10~ 2 

< 1.0 x 10~ 2 

< 4.4 x 10~ 3 

< 5.8 x 10~ 2 

6.0 x 10- 3 

< 8.2 x 10' 2 

so/c 18 o 

1.8 x 1(T 3 

3.3 x 10~ 3 

3.0 x 10- 3 

2.1 x 10~ 3 

1.7 x 10~ 3 

< 1.1 x 10~ 2 

< 7.1 x 10~ 4 

< 1.0 x 10' 2 

c-C 3 H 2 /C 18 0 

< 3.4 x KT 1 

2.0 x 10~ 3 

5.1 x 10- 3 

< 3.3 x 10~‘ 

< 2.1 x 10- 1 

< 6.1 x 10“* 

< 7.1 x 10- 2 

< 1.1 x 10- 2 

hco/c 18 o 

< 4.5 x 1(T 4 

1.2 x 10~ 3 

2.0 x 10“ 3 

< 1.3 x 10~ 3 

2.0 x 10~ 3 

< 1.1 x 10~ 2 

< 9.2 x 10~ 4 

< 1.6 x 10~ 2 

c 2 h/c 18 o 

4.5 x 1(T 2 

5.5 x 10~ 2 

1.4 x 10“* 

3.1 x 10~ 2 

8.5 x 10- 2 

1.6 x 10- 1 

3.5 x 10- 2 

6.9 x 10~ 2 

n 2 h + / c 18 o 

5.9 x 1(T 4 

7.8 x 10~ 4 

5.1 x 10~ 4 

7.9 x 10- 4 

2.5 x 10- 3 

< 1.3 x 10~ 3 

9.5 x 10- 4 

< 1.1 x 10~ 3 

hnc/c 18 o 

2.8 x 10~ 3 

3.5 x 10~ 3 

4.4 x 10“ 3 

2.6 x 10~ 3 

5.4 x 10~ 3 

6.6 x 10- 3 

5.2 x 10“ 3 

5.1 x 10~ 3 

hcn/c 18 o 

7.8 x 1(T 3 

7.5 x 10' 3 

1.1 x 10“ 2 

8.6 x 10- 3 

1.9 x 10~ 2 

1.6 x 10~ 2 

7.4 x 10~ 3 

1.7 x 10~ 2 

CN / C 18 0 

2.9 x 1(T 2 

3.4 x 10~ 2 

3.9 x 10- 2 

2.5 x 10- 2 

5.3 x 10- 2 

7.2 x 10~ 2 

9.8 x 10~ 3 

3.7 x 10- 2 

HCN / HCO + 

1.9 

2.0 

1.1 

3.4 

2.7 

1.4 

4.1 

2.6 

CN / HCN 

3.8 

4.5 

3.5 

2.9 

2.8 

4.4 

1.3 

2.1 

HCN/CS 

1.5 

0.9 

1.2 

2.5 

3.2 

1.3 

1.4 

2.4 


The errors of the ratios range in general between 2% and 63% and have an average value of 14%. Exceptions to these values are the ratios 
involving HNCO and CH 3 CN in M 82, as their column densities have large errors associated. 


telescope. Thanks to the smaller beam size of that telescope 
at A -3 mm, their observations can separate the emission of the 
CND and that of the starburst ring, and show that the emission of 
some species such as C 2 H and N 2 H + come mainly from the lat¬ 
ter. Unfortunately, there are no similar estimations for NGC 7469 
and M 51 calculated for the IRAM 30 m beam at the 3 mm wave¬ 
length of our survey. 

CH 3 CCH and C-C 3 H 2 were not detected in any of the AGNs, 
while HCO and HOC + are only tentatively seen in NGC 1068, 
and have similar abundances than in the starburst galaxies. 
Regarding shocks and dust grain chemistry, most of the tracers 
(such as OCS, HNCO, or SO) are only detected in NGC 1068. 
Shocks seem to be affecting part of the ga s in the CND of this 
galaxy (iKrips et al.ll201 It IViti et al.ll2014l) . CH 3 OH is also de¬ 
tected in M 51 (not in NGC 7469), but with a clearly lower abun¬ 
dance than in NGC 1068. 


Among others, we examined the following relative abun¬ 
dances: HCN/C ls O, CN/C ls O, HCN/HCO + , CN/HCN, and 
HCN/CS (Table 0). The last three were proposed as dis¬ 
criminators between AGNs and starburst activities based 


on two-dimens i onal p l ots of line ratios ( iKohno et al] [2001: 
Meii erink et alJ 120071 : iKrins et all 120081 : llzumi et ak l2013h . 


Using our data, HCN/CS, HCO+/HCN, and CN/HCN, either by 
themselves and plotting one ratio versus other, do not show a 
correlation with the galactic nuclear activity of our sample. The 
most probable reason for the lack of trend between all these line 
ratios and AGNs (versus starburst galaxies) is likely the low an¬ 
gular resolution of the IRAM 30 m telescope, which is not suffi¬ 
cient to separate the CND and the surrounding starburst regions. 
However, we do see that the HCN abundance is clearly enhanced 
in AGNs, and the reason for that_could be a high temperature- 
driven chemistry (llzumi et al.ll2013l) . or shocked regions at a few 
hundred parsecs from the super mas sive black holes due to out¬ 
flowing material (Martin et all 2015). Interestingly, HCN/C ls O, 
CN/C ls O, C 2 H/C 18 0, and HNC/C ls O in M51 are closer to the 
starburst results than to the other AGNs (Fig. [ 6 }. 


SiO is tentatively detected in NGC 1068 but not in M51 
and NGC 7469. In the case of NGC 7469, where lines are 
weaker, we cannot conclude from the upper limit to the SiO 
abundance whether the gas chemistry is different. This species 


usually traces shocks in Galactic and extragalactic sources 
(Martin-Pintad o et al. _ 1992 19971: Garcla-Burillo et all [2000, 


12001 : Riquelme et al. 20 lot Takano et al.l2014l) . However, some 

studies (based on observational data or on chemical models) 
point to SiO as a tracer of X-rays an d cosmic rays in the centres 
of the Milky Way and NGC 1068 (Martin-Pintado et afl [2000; 


lUsero et alJl2004l : lAladro et alJl2013h . If that is the case in all 

AGNs, its detection in the other two Seyferts of our sample 
would be somehow expected (if SiO was strong enough to be 
detectable). A study devoted to this molecule is necessary to bet¬ 
ter understand how it is formed under the different regimes (i.e. 
shock-dominated versus X-ray dominated gas). 


Using the IRAM 30 m telescone. lWatanabe et al.l (120141) car¬ 
ried out a molecular line survey, in the 3 mm and 2 mm bands, 
outside the nucleus of M 51. They observed two regions in a 
spiral arm of the galaxy. Their position PI has twice the SFR 
of position P2. Both SFRs are quite low, in the range [0.02- 
O.O5]M 0 yr _1 . This allows us, for the first time, to compare the 
nuclear molecular gas composition with that outside the central 
region. Comparing the fractional abundances (with respect to 
C 18 0) of the twenty-two species in common between our work 
and Watanabe’s, we find that the chemical composition of the 
gas is fairly the same, except for a few cases. H 13 CO + , HNCO, 
and C-C3H2 are detected in the PI position, but not in the cen¬ 
tre. On the other hand, we detect H I 3 CN, C 34 S, and SO, which 
are not seen in the P2 position. Apart from that, the rest of the 
molecules have similar abundances, within a factor of <2 in the 
three positions. One would expect that the chemistry of the gas 
might change between the nucleus, affected by higher star for¬ 
mation and the proximity of the black hole, and the spiral arms. 
However, given the beam size of the IRAM 30 m at these fre¬ 
quencies and the source sizes of the three positions, the emis¬ 
sion might be so diluted that the differences are cancelled out. 
Further investigations with interferometers are needed to bet¬ 
ter understand the change of the gas chemistry as a function of 
galactocentric radius. 
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7.3. ULIRGs 

Independently of their powering source, the ULIRGs of our sam¬ 
ple show some chemical differences with respect to the AGNs 
and starburst galaxies described in the previous sections, as well 
as some dissimilarities between Arp 220 and Mrk231 them¬ 
selves. 

The HNC (1-0) line intensity in Arp 220 slightly exceeds 
that of HCN(1 - 0) by a factor of 1.3, while for Mrk231 we 
obtain HNC/HCN (1 - 0) — 0. 6. This particular high rat io in 
Arp 220 was alrea dy noticed by Huettemei ster et al.l (1 19951) and 
Aalto et al.l (1 20071) . and seems to be due to an overabundance of 
HNC in this galaxy, a consequence of either mid-IR pumping 
of the lines, or a chemistry driven by X- ray dominated regions 
(Meiie rink et al.l [20071 lAalto et aTI 2007). A priori, we cannot 
rule out any of the two scenarios with our data, since we would 
need several HCN and HNC transitions in order to estimate the 
H 2 column and volume densities, indispensable to address the 
situation. However, we dare to favour IR pumping effects, given 
that the chemical composition of Arp 220 is not similar to the 
AGNs of our sample (XDR scenario), as we describe below. 

CH 3 CN and SiO are relatively abundant in Arp 220 
(fevvxlO 5 with respect to C ls O), but are not detected in 
Mrk 231. On the other hand, HCO + , H°CN, and CN are four 
times more abundant in Mrk 231 than in Arp 220. Our results 
suggest that the Arp 220 chemistry resembles more the starburst 
galaxies of our sample dominated by shocks, with high abun¬ 
dances of HC 3 N, CH 3 OH, HNCO, CH 3 CN and SiO relative to 
C ls O (10~ 2 ~ 3 ), but no PDR tracer detections (C-C 3 H 2 , HCO, 
HOC + ), except for CH 3 CCH (just tentatively detected). This 
would support the infrared pumping of HNC in the galaxy at the 
expense of a XDR environment, as mentioned before. Mrk231, 
in contrast, is characterised by high abundances of CN, HCN and 
H 13 CN relative to C 18 0 (10 2 ' 3 ), more similar to those found in 
NGC1068 and NGC 7469. 

Vibrationally excited lines from HC 3 N, arising from dense 
(n > 10 6 cm~ 3 ) and hot (T ex = 190 + 20 K) cores, are detected 
in Arp 220. They are a factor 4-5 fainter than the correspond¬ 
ing rotational transitions. Assuming that this ratio may be simi¬ 
lar for the rest of the galaxies, the expected vibrational intensi¬ 
ties in Mrk 231 would be ~0.2 mK, which is well below our rms 
level. Thus, we cannot discard lower level vibrational emission 
in Mrk 231. However, vibrationally excited HC 3 N lines would 
be observable in NGC 253 and M82, but are not detected. This 
clearly points to a different excitation between normal starbursts 
and ULIRGs. 

There are also common features between the Mrk 231 and 
Arp 220 spectra as compared t o the rest of the sample. In 
agreement with previous s tudies dGreve et all200 9: Martin et al J 
1201 It iHenkel et al.l 120141) 13 CO and C 18 0 are equally abun¬ 
dant in both ULIRGs. This may be due to high opacities of 
the 13 CO lines, to an intrinsic underabundance of l3 CO in lumi¬ 
nous infrared galaxies (Tani guchi et al.l[l999l) . or/and to a C 18 Q 
overabundance due to t he strong starbursts jMatsushit a et alJ 
l2009t iMartfn et al.ll201 ll) . High opacities of 13 CO in these two 
galaxies seem unlikely, given the weakness of the 13 CO(l - 
0) lines relative to 12 CO(l - 0) (they are >20 times fainter). 
An overabundance of l8 0 i n Arp 220 has been recently found 
by iGonzalez-Alfonso et al.l ( 2014 ). However, the 12 C/ 13 C and 
16 0/ 18 0 line ratio s obtained for Arp 220 and Mrk231 by 
iHenkel et al.l ( 2014 ) do not favour this last scenario. It can be 
seen in Fig. [ 6 ]that the 13 CO/C ls O ratios in Mrk231 and Arp 220 
are a factor of ~4 lower than in the rest of the galaxies. We ob¬ 
tain a 12 CO/ 13 CO(l - 0) relative intensity of 21+7, which is a 


factor two lower than that obtained by iGreve et ail ( 2009 ) for 
Arp 200 (43+10). We note that our 12 CO (1 -0) and 13 CO (1 -0) 
lines were observed simultaneously with the EMIR receiver of 
the 30m telescope, while IGreve et all ( 2009 ) used older data to 
derive this ratio, which may be less accurate. 

In addition, HC 3 N/C 18 0 ratios are clearly higher in the 
ULIRGs than in the more nearby starbursts and AGNs ob¬ 
served by us. If C 18 0 is overabundant in ULIRGs, then HC 3 N 
should be it as well. This effect was already observed in other 
LIRGs/ULIRGs, and the reason could be a higher amount of 
dense gas and warm d ust that protects the species from be¬ 
ing photo-disso ciated (lAalto et al. 1120071: ICostagliola et alJl 201 lb 
Lindb erg et all 1201 ll) . Thus, HC 3 N seems to be a well suited 
species to study the activity in ULIRGs, regardless of the power¬ 
ing source in their nuclei. Also common to Arp 220 and Mrk 231 
is the outstanding H 13 CN/C 18 0 ratio, as previously commented 
in Sect. [ 6 ] Again, if C 18 0 is overabundant in Arp 220 and 
Mrk 231, it only makes more clear the difference in abundances 
of HC 3 N and H 13 CN between the ULIRGs and the rest of the 
sample. On the other hand, if C ls O has similar values in all the 
galaxies, it results in an even bigger contrast. 

It is also interesting to compare ULIRGs, characterised by 
strong interactions between galaxies, with isolated galaxies, 
where the chemical and physical properties are the result of in¬ 
trinsic and secular mechanisms. In CIG638, a_good representa¬ 
tive of an isolated galaxy (see IVerlev et al.l (20071 for the crite¬ 
ria used to define these objects), the C 2 H/HCN and C 2 H/HNC 
integrated intensity ratios were found to be interestingly high 
(1.0 and 4.3 respectively), either because the galaxy has a 
lesser amount of dense gas, or/and due to an overproduction 
of C 2 H (iMartln et al.ll2014l) . The C 2 H/HCN ratios for Arp 220 
and Mrk 231 are a factor of 3 and 4 (respectively) lower than in 
CIG638. Moreover, the C 2 H/HNC integrated intensity ratio is 
even more contrasted, being 9 and 5 times (respectively) lower 
in the ULIRGs than in the isolated galaxy. This points to a very 
likely difference in the gas composition between these two types 
of galaxies, with the ULIRGs containing more dense gas (asso¬ 
ciated with their nuclear starbursts), than the isolated galaxies. 
On the other hand, the rest of the galaxies of our sample have a 
C 2 H/HNC ratio in the range [0.7 -1.5] which, according to the 
results derived by Martin et al] d2014l) . agree with the expected 
values for “normal” (local) starbursts and AGNs. These results 
are on average well below the CIG638 and above the Arp 220 
and Mrk231 values. 

8. Summary and conclusions 

Using the IRAM 30 m telescope, we performed a molecular 
line survey towards the circumnuclear regions of eight active 
galaxies, covering the frequency range —[86-116] GHz. The 
sample is composed by three starburst galaxies, M83, M82, 
and NGC 253, three AGN galaxies with starburst influence, 
NGC 1068, M51, and NGC7469, and two ULIRGs, Mrk231 
and Arp 220. The column densities of twenty-seven species (or 
their upper limits when they are not detected) were calculated 
assuming LTE and a single source size for each galaxy based 
on previous available interferometric maps. The studied species 
include ten 13 C, 18 0, 17 0, and 34 S bearing isotopologues. We 
compared the abundances, with respect to C 18 0, among the eight 
galaxies to characterise the average molecular gas composition 
of each kind of galaxies. Our results are summarised in Figs. Q] 
and [ 6 ] and Tables [2] 0 and [4] 

The comparison among the galaxies allowed us to find out 
the key species that show the most contrasted abundances, rel- 
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ative to C ls O, among the sources. CH 3 CCH is only detected in 
starbursts and is a factor of one order of magnitude more abun¬ 
dant in M82 (the galaxy with the highest abundance) than in 
NGC 1068 (not detected). This molecule seems to be associated 
with massive star formation, and might be a good tracer of PDRs. 
H 13 CN, CH 3 CN and SiO have a peak of abundance in Arp 220, 
which is twenty-one (H 13 CN), and seven (CH 3 CN, SiO) times 
respectively higher than in M 82 (minimum abundance). We find 
high C 2 H abundances in M 82, NGC 1068 and NGC 7469, being 
in NGC 7469 five times higher than in M 51 (minimum). A large 
part of the C 2 H emission in the AGNs is expected to come from 
the starburst rings that surround the CNDs. HNCO shows abun¬ 
dance variations among the starburst galaxies of a factor of four. 
This latter result supports previous claims of HNCO a s a g ood 
tracer of the starburst evolutionary state (iMartfn et alJl2009al) . 
Other species, such as HCN and HC 3 N, show smaller, though 
still significant, abundance differences among galaxies by fac¬ 
tors of three to four. All the cited molecules constitute good can¬ 
didates to study the chemistry in galaxy nuclei. 

Regarding the starburst galaxies, the molecular gas of M83 
is characterised by relatively high abundances (with respect 
to C 18 0, 10~ 2 -1(U 4 ) of some shock tracers, such as CH 3 OH, 
HNCO and SiO, as well as the absence of both Ha recombi¬ 
nation lines and the PDR tracers C-C 3 H 1 and HCO. In NGC 253 
shocks, created by molecular cloud collisions and stellar winds, 
dominate the heating of the gas, as shown by the presence 
of several other shock- and dust grain-tracers (e.g. SO, S0 2 , 
OCS, HOCO + ). In addition, the molecular cloud complexes are 
ionised by OB clusters, which translates into high intensities of 
Her recombination lines, as well as moderate abundances of c- 
C 3 H 2 , HOC + , and HCO (KT 3 -1(T 4 ). Lastly, M82 represents a 
galaxy where UV fields are very intense and create large PDRs. 
Her recombination lines, as well as C-C 3 H 1 , HCO, HOC + and 
CH 3 CCH are brighter than in the other starbursts of the sample. 
CH 3 OH and HNCO are still observable in the cores of the molec¬ 
ular clouds, although their abundances are clearly lower than in 
NGC 253 and M 83, due to the lack of shielding from UV pho¬ 
tons. We explore a possible scenario where M 83, NGC 253 and 
M 82 are considered templates of galaxies in young, intermedi¬ 
ate and late phases of starbursts, and follow the evolution of the 
gas composition along time. 

As for the AGNs (NGC 1068, NGC 7469, M51 and 
Mrk 231), HCN/HCO+ and HCN/CS are quite low in NGC 7469, 
similar to the values in starburst galaxies; HCN/C ls O and 
CN/C ls O, C 2 H/C 1 s O, and HNC/C ls O are also much the same 
in starbursts and in the low luminosity AGN of M 51. The lack 
of molecular trends in our AGN galaxies, as compared to the 
starbursts, seems to be due to their CNDs not being spatially re¬ 
solved by the IRAM 30 m beam. Remarkably, CH 3 OH, SiO and 
HNCO are only detected in NGC 1068. 

Finally, the ULIRGs Mrk 231 and Arp 220 show common 
spectral features, such as the almost identical 13 CO and C ls O 
column densities, and higher H 13 CN/C 18 0 and HC 3 N/C ls O 
abundance ratios than the rest of the galaxies. Nevertheless, both 
ULIRGs differ in some line detections that seem to point out to 
different mechanisms at work in their nuclear regions. In particu¬ 
lar, Arp 220 shows emission of HC 3 N vibrationally excited lines 
that arise from warm gas (7\.j h = 190 + 20 K), and shock tracers 
such as CH 3 OH, HNCO, and SiO. This, together with the lack of 
PDR tracers such as HOC + , HCO, and c-C 3 H 2 , makes Arp 220 
resemble the spectrum of a giant hot core-like starburst. On the 
other hand, Mrk231 has higher abundances of HCN, H 13 CN, 
and HCO + , like the values of NGC 1068 and NGC 7469, thus 
being more similar to an AGN-dominated galaxy. 
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Fig. 1: Spectra of all eight galaxies between 86 and 116 GHz. The strongest molecules are indicated. The faint detected species are 
shown in Figs.|2]and [3] The spectra refer to the rest frame adopting the velocities in Table[T| 
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Fig. 2: Zoomed spectra between 86 GHz and 100 GHz where the faint lines detected are highlighted. Lines in brackets were not 
detected in all the galaxies (typically they were only detected in NGC253). For details about the detections see Appendix [A] The 
spectra refer to the rest frame adopting the velocities in Table Q] 
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Aladro et al. (2015): A =3 mm line survey of nearby active galaxies 
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Fig. 3: Continued zoomed spectra between 100GHz and 116 GHz where the faint lines detected are highlighted. Lines in brackets 
were not detected in all the galaxies (typically they were only detected in NGC 253). For details about the detections see AppendixlAl 
The spectra refer to the rest frame adopting the velocities in Table [T] 
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Fig. 4: Lowest frequency range observed exclusively towards NGC253 and M82 (the rest of the galaxies were observed from 
~ 86 GHz). 
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Fig. 5: Zoomed spectra of NGC 253 showing the Gaussian fits to the three strongest transitions of HC 5 N detected tentatively. Dashed 
lines are lcr and 3 <x noise levels. 
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Fig. 6: Abundances of thirty-four species relative to C 18 0 for the eight galaxies of our sample. Dashed histograms indicate upper 
limits. 
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Table A. 1: Results for M 83 


Line 

Frequency 

MHz 

/ 

Kkms~' 

ff.SR 

km s _1 

AT 

km s _1 

■^Mtf Comments 
mK 

H 13 CN (1 -0) 

86340.2 

0.25 ± 0.04 

503.4 ± 6.8 

64.7 ± 16.0 

3.55 

hf. t 

H l3 CO + (l -0) 

86754.3 

0.21 ±0.06 

510 

95 

2.07 

t 

SiO (2-1) 

86847.0 

0.24 ± 0.06 

510 

95 

2.40 

t 

HN I3 C(1 -0) 

87090.8 

0.11 ±0.05 

510 

95 

1.11 

t 

C 2 H(1 -0) 

87316.9 

4.01 ± 0.08 

504.3 ± 1.3 

91.7 ±3.2 

22.62 

hf 

HNCO (4o.4 - 3 0> 3) 

87925.2 

0.69 ± 0.09 

510 

95 

6.79 

m 

HCN (1-0) 

88631.8 

10.18 ±0.13 

501.7 ±0.8 

95.5 ± 1.8 

97.47 

hf 

HCO + (1-0) 

89188.6 

9.30 ±0.15 

502.8 ± 1.0 

95.0 ± 2.4 

90.96 


HOC + (1 - 0) 

89487.4 

0.29 ± 0.04 

516.3 ±7.7 

95 

2.87 

t 

HNC (1-0) 

90663.6 

4.12 ±0.08 

505.3 ± 1.3 

93.0 ±3.0 

41.45 


HC 3 N(10-9) 

90979.0 

0.44 ± 0.05 

512.8 ±6.4 

95 

4.34 


CH 3 CN(5 a - 4 a ) 

91987.0 

0.24 ± 0.06 

510 

95 

2.33 

m, t 

N 2 H + (1 - 0) 

93173.7 

1.17 ±0.05 

510 

95 

13.32 

m 

C 34 S (2-1) 

96412.9 

0.42 ± 0.05 

510 

95 

4.19 


CH 3 OH(2 A -l A ) 

96741.4 

1.08 ±0.05 

504.1 ± 1.9 

57.5 ±4.5 

17.13 

m 

CS (2 - 1) 

97981.0 

3.25 ± 0.07 

501.9 ± 1.4 

95 

32.04 


SO (2, - 3 2 ) 

99299.9 

0.45 ± 0.05 

503.3 ± 9.3 

120.0 ±22.1 

3.49 

t 

HC 3 N(11 - 10) 

100076.4 

0.36 ± 0.04 

512.8 ±6.4 

95 

3.51 

t 

CH 3 CCH(6 A -5 A ) 

102548.0 

0.44 ± 0.05 

521.5 ±22.2 

95 

4.23 

m, t 

CH 3 OH(0 0 - l-i) 

108893.9 

0.20 ± 0.01 

504.1 ± 1.9 

57.5 ±4.5 

3.28 

t 

HC 3 N(12- 11) 

109173.6 

0.26 ± 0.03 

512.8 ±6.4 

95 

2.61 

t 

C 18 0(1 -0) 

109782.2 

3.40 ± 0.07 

502.2 ± 1.3 

94.7 ±3.1 

33.65 


HNCO(5 0 .5-4o, 4 ) 

109905.8 

0.70 ± 0.09 

510 

95 

6.87 

m 

13 CO(l -0) 

110201.4 

17.29 ± 0.24 

502.2 ± 0.9 

97.8 ±2.1 

165.00 


CH 3 CN(6 A -5 A ) 

110383.5 

0.27 ± 0.06 

510 

95 

2.65 

m, t 

C 17 0(1 -0) 

112359.3 

0.63 ± 0.06 

510 

95 

6.24 

t 

CN (l 0 j - 0 0 ,i) 

113191.3 

6.31 ±0.06 

500.7 ± 0.9 

118.4 ±2.2 

31.14 

m 

CN (1 0 , 2 - 0 0 ,i) 

113491.0 

12.66 ±0.13 

500.7 ± 0.9 

118.4 ±2.2 

91.81 

m 

12 CO(l -0) 

115271.2 

225.92 ± 2.85 

504.4 ± 1.3 

79.6 ± 2.8 

1940.18 



Remarks: ( b ) blended line; (m) multi-transition line; (hf) hyperfine transition; (f) tentative detection. Parameters without errors were fixed 
when fitting the lines. 


Appendix A: Simulated fit results for individual spectra. 


Table A.2. Results for NGC 253 


Line 

Frequency 

MHz 

/ 

Kkms - 1 

Flsr 

kms~' 

AV 
km s ' 1 

r^-peak 

1 MB 

mK 

Comments 

NH 2 CN(4 u -3,, 2 ) 

80504.6 

0.70 + 0.15 

250 

180 

3.68 

t 

C 2 S ( 6 7 - 5 6 ) 

81505.2 

0.66 + 0.18 

250 

200 

3.08 

t 

HC 3 N (9 - 8 ) 

81881.5 

7.37 + 0.15 

250 

180 

38.41 


c-C 3 H 2 (2 0 , 2 - 1 1 , 1 ) 

82093.6 

1.66 + 0.07 

250 

180 

8.67 


HC 5 N (31-30) 

82539.3 

0.99 ± 0.45 

240 

180 

5.15 

t 

c-C 3 H 2 (3i j2 - 3 o, 3 ) 

82966.2 

0.76 ± 0.03 

250 

180 

3.94 

t 

CH 3 OH (5_i - 4o) 

84521.2 

2.52 + 0.12 

249.4 + 5.1 

217.7+11.9 

10.9 


OCS (7 - 6 ) 

85139.1 

1.39 + 0.08 

250 

200 

6.53 

b 

HC ls O + (1 - 2) 

85162.2 

0.84 ± 0.04 

250 

200 

3.93 

b 

HC 5 N (32-31) 

85201.6 

0.92 ± 0.41 

240 

180 

4.78 

b 

c-C 3 H 2 (2i j2 - l 0 ,i) 

85338.9 

5.99 ± 0.26 

250 

180 

31.09 


CH 3 CCH (5k - 4a;) 

85457.3 

2.66 ± 0.29 

250 

200 

12.39 

m 

HOCO + ( 4 0 ,4 - 3o, 3 ) 

85531.5 

0.38 ± 0.22 

250 

180 

1.98 

t 

H42a 

85688.4 

- 

230 

180 

8 

b 

SO(2 2 -10 

86094.0 

0.88 ± 0.07 

250 

180 

4.60 


H 13 CN(1 -0) 

86340.2 

4.71+0.13 

250.1 +3.2 

181.5 + 7.4 

24.25 

hf 

HCO(li,o-Oi, 0 ) 

86670.8 

1.92+ 1.17 

250 

200 

4.68 

b, hf 

H 13 CO + (1 - 0) 

86754.3 

3.81+0.44 

250 

200 

17.86 

b 

SiO (2- 1) 

86847.0 

4.76 ± 0.44 

250 

200 

22.30 


HN 13 C(1 -0) 

87090.8 

2.09 + 0.12 

250 

200 

9.81 



Continue in next page 
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Continued from previous page 
C 2 H(1 -0) 

HC 5 N (33-32) 
HNCO(4o,4-3o, 3 ) 

HCN (1-0) 

HCO + (1-0) 

HOC + (l - 0) 

HNC (1 - 0) 

HC 5 N (34-33) 

HC 3 N (10 - 9) 

CH 3 CN (5* - 
H41a 

13 CS (2 - 1) 

N 2 H + (1 - 0) 

HC 5 N (35-34) 

CH,CHO(5i 5 -4 i4 )A 
CH 3 CHO(5i, 5 -4 1i4 )E 
C 2 S (7 g - 6 7 ) 

CH 3 OH (80-70 
HC 5 N (36-35) 

CH 3 OH(2i - li) 
CH 3 CHO(5 o,5-4o, 4 )E 
CH 3 CHO(5o, 5 -4o, 4 )A 
C 34 S (2-1) 

CH 3 OH(2,-U) 

OCS (8 - 7) 

CH 3 OH(2, - li) 

CS (2 - 1) 

HC 5 N (37-36) 

H40q' 

SO (2 3 - 1 2 ) 

NH 2 CN(5i, 5 -4 m ) 

NH 2 CN(5o,5-4 0 , 4 ) 

HC 3 N(11 - 10) 
NH 2 CN(5 w -4 1>3 ) 

HC 5 N(38 - 37) 

H 2 CS (3i, 3 -2 u ) 
CH 3 CCH(6,-5,) 

H 2 CS (3o, 3 - 2 0 , 2 ) 

HC 5 N(39- 38) 

SO 2 (3i i3 - 2o, 2 ) 

H 2 CS (3i, 2 -2 u ) 

C 2 S (89 - 7g) 

H39fr 

HOCO + (5o,5-4 0 , 4 ) 

13 CN(1 2j i - 0i, 0 ) 
13 CN(1 u - 0 u ) 
13 CN(1 2 ,2-0 u ) 

CH 3 OH(0 0 - l-i) 

HC 3 N (12 - 11) 

SO(3 2 -20 
OCS (9 - 8 ) 

C ls O(l -0) 
HNCO(5o,5-4 0 , 4 ) 

13 CO(l -0) 

CH 3 CN(6t-5*) 

C n O(l -0) 

CN(l 0 ,i-Oo,i) 

CN(1o, 2 -Oo,i) 

NS(3i,_i- 2 u ) 

12 CO(l -0) 

H38tr 

Continue in next page 
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87316.9 

34.64 ± 0.45 

246.9+ 1.9 

186.8 + 3.8 

106.10 

hf 

87863.9 

0.84 + 0.38 

240 

180 

4.39 

t, b 

87925.2 

5.79 + 1.73 

250 

200 

26.99 

m, b 

88631.8 

68.37 + 0.72 

249.1 ± 1.3 

199.7 + 3.1 

305.20 

hf 

89188.6 

57.92 + 0.63 

254.8 + 1.4 

199.7 + 3.4 

266.77 


89487.4 

2.13+0.11 

250 

200 

10.02 


90663.6 

35.45 ± 0.27 

238.6 + 0.9 

185.1+2.2 

177.46 

b 

90526.2 

0.76 ± 0.34 

240 

180 

3.99 

t, b 

90979.0 

8.42 + 0.18 

250 

180 

43.90 


91987.0 

3.59 + 2.66 

250 

180 

18.55 

m, b 

92034.4 

- 

230 

180 

5 

b 

92494.3 

0.99 ± 0.07 

250 

200 

4.65 


93173.7 

10.64 + 0.19 

235.7 + 2.0 

180 

54.89 

m, b 

93188.5 

0.69 + 0.31 

240 

180 

3.59 

t, b 

93580.9 

0.52 ± 4.05 

250 

200 

2.44 

t 

93595.2 

0.52 ± 4.03 

250 

200 

2.43 

t 

93870.1 

0.68 + 0.19 

250 

200 

3.18 

t 

95169.4 

0.98 + 0.12 

233.3 + 10.6 

171.8 + 24.2 

5.38 


95850.7 

0.61+0.28 

240 

180 

3.21 

t, b 

95914.3 

0.81+0.01 

250 

200 

3.80 

b 

95947.4 

0.65 ± 5.03 

250 

200 

3.04 

b, t 

95963.4 

0.65 ± 5.06 

250 

200 

3.06 

b, t 

96412.9 

4.39 + 0.13 

250 

200 

20.58 


96741.4 

17.14 + 0.26 

250 

200 

78.40 

m, b 

97301.2 

1.38 + 0.08 

250 

200 

6.47 


97582.8 

0.83 + 0.01 

250 

200 

3.89 


97981.0 

33.98 ± 0.46 

234.5+ 1.7 

195.2 + 4.1 

162 


98512.9 

0.54 ± 0.24 

240 

180 

2.83 

t 

99023.0 

1.45 + 0.50 

237 

180 

7 


99299.9 

4.20 ± 0.32 

250 

180 

21.91 

b 

99311.2 

1.07 + 0.23 

250 

180 

5.58 

m, b 

99972.7 

1.06 + 0.23 

250 

180 

5.51 


100076.4 

9.02 + 0.19 

250 

180 

47.02 


100629.5 

1.09 + 0.24 

250 

180 

5.66 

m 

101175.1 

0.48 + 0.21 

240 

180 

2.49 

t 

101477.6 

0.58 + 0.12 

250 

180 

3.02 

t 

102548.0 

3.89 ± 0.42 

250 

200 

18.10 

m 

103040.2 

0.83 + 0.18 

250 

180 

4.33 

t 

103837.3 

0.41+0.19 

240 

180 

2.16 

t 

104029.4 

0.56 ± 0.32 

240 

180 

2.90 

t 

104617.0 

0.60 + 0.13 

250 

180 

3.15 

t 

106349.7 

0.61+0.17 

250 

200 

2.87 

t 

106737.4 

1.75 + 0.34 

220 

181 

9 


106913.6 

0.59 ± 0.34 

250 

180 

3.07 

t 

108651.3 

1.08 + 0.09 

250 

200 

9.72 

m, b 

108657.6 

1.08 + 0.09 

250 

200 

9.72 

m, b 

108780.2 

1.94 + 0.16 

250 

200 

8.92 

m 

108893.9 

3.50 ± 0.05 

250 

200 

16.39 


109173.6 

9.12 + 0.19 

250 

180 

47.56 


109252.2 

1.40 + 0.11 

250 

180 

7.28 

t 

109463.1 

1.20 + 0.07 

250 

200 

5.65 

b, t 

109782.2 

22.49 ± 0.57 

258.3 + 3.3 

204.5 ± 7.8 

102.65 


109905.8 

7.62 ± 2.28 

250 

200 

35.55 

m 

110201.4 

81.16 + 0.96 

259.9+ 1.6 

200.8 ± 3.6 

370.67 


110383.5 

5.03 ± 3.73 

250 

180 

25.94 

m 

112359.3 

3.15 + 0.14 

269.9 + 5.7 

200 

14.77 


113191.3 

47.93 ± 0.32 

249.6 ± 0.8 

193.2 + 2.0 

184.30 

m 

113491.0 

95.53 ± 0.62 

249.6 ± 0.8 

193.2 + 2.0 

425.50 

m 

115153.9 

1.80 + 0.90 

250 

180 

8.86 

m, b, t 

115271.2 

1069.30+8.60 

260.9 + 0.1 

210.6 + 0.1 

4769.30 

b 

115274.4 

- 

230 

180 

7 

b 
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Continued from previous page 


NS(3 U -2 1 ,- 1 ) _ 

See caption of Table lA. ll for details. 


115556.2 


1.80 + 0.90 


250 


180 8.32 m, t 


Table A.3: Results for M 82 


Line 

Frequency 

MHz 

/ 

K km s 1 

Vlsr 

kms -1 

AV 
km s 1 

7C Comments 
mK 

C-C 3 H 2 ( 32 , 2 - 3 !, 3 ) 

84727.7 

0.143 + 0.06 

300 

100 

1.34 

t 

HC ls O + (l -0) 

85162.2 

0.41 ±0.06 

300 

100 

3.82 

b 

C-C 3 H 2 ( 2 12 - 1 1>0 ) 

85338.9 

3.36 + 0.13 

300 

100 

31.16 


CH 3 CCH(5 k -4 k ) 

85457.3 

2.59 ± 0.08 

314.6+1.9 

100 

23.94 

m 

H42a 

85688.4 

- 

310 

100 

11 

b 

H 13 CN(1 -0) 

86340.2 

0.38 ± 0.04 

300 

100 

3.54 

hf, t 

HCO di , 0 - 0i, 0 ) 

86670.8 

0.32 ± 0.03 

300 

100 

2.99 

b, t 

H 13 CO + (1-0) 

86754.3 

0.89 ± 0.04 

300 

100 

8.39 

b 

SiO(2 - 1) 

86847.0 

0.30 ± 0.04 

300 

100 

2.87 

b, t 

C 2 H(1 -0) 

87316.9 

20.61+0.23 

294.1 ± 1.1 

134.2 ± 2.4 

83.31 

hf 

HNCO( 4 0 , 4 - 3 o, 3 ) 

87925.2 

0.26 ± 0.45 

300 

100 

2.41 

m, t 

HCN (1 - 0) 

88631.8 

23.56 ± 0.40 

294.1 ± 1.1 

113.6 + 2.6 

182.00 

hf 

HCO + (1-0) 

89188.6 

33.74 ± 0.80 

292.8 ± 1.6 

128.8 + 3.7 

236 


HOC + (1 - 0) 

89487.4 

0.56 ± 0.04 

300 

100 

5.21 


HNC(l-O) 

90663.6 

10.48 ± 0.23 

296.3 ± 1.6 

115.8 + 3.8 

83.79 


HC 3 N (10-9) 

90979.0 

1.42 + 0.03 

305.7 ± 1.5 

109.5 ± 3.4 

11.64 


CH 3 CN(5 t -4*) 

91987.0 

0.24 ± 0.72 

300 

100 

2.19 

m, b 

H41a 

92034.4 

- 

280 

120 

9 

b 

N 2 H + (1 -0) 

93173.7 

1.65+0.04 

299.5 ± 1.5 

100 

15.35 

m 

CH 3 CHO (5i j5 - 4 i 4 ) A 

93580.9 

0.26 ± 0.02 

300 

100 

2.12 

m, t 

C 34 S (2-1) 

96412.9 

0.58 ± 0.03 

294.9 ± 3.3 

100 

5.46 


CH 3 OH (2 k — 1*) 

96741.4 

1.01 ±0.03 

300 

100 

9.17 

m 

CS (2 - 1) 

97981.0 

9.97 + 0.16 

294.8 ±1.1 

109.2 ± 2.6 

84.66 


H40tr 

99023.0 

1.40 + 0.70 

302 + 8 

100 

13 


SO (2, -3 2 ) 

99299.9 

1.11 ± 0.11 

300 

100 

10.39 


HC 3 N(11 - 10) 

100076.4 

1.42 ±0.03 

305.7 + 1.5 

109.5 ± 3.4 

12.26 


CH 3 CCH (6k - 5t) 

102548.0 

3.48 ±0.11 

314.6 ± 1.9 

100 

32.01 

m 

H39a 

106737.4 

1.00 + 0.50 

307.2 ± 5.4 

100 

9 


CH 3 OH(0 0 - l_i) 

108893.9 

0.254 ± 0.009 

300 

100 

2.39 

m 

HC 3 N(12- 11) 

109173.6 

1.34 ± 0.03 

305.7 ± 1.5 

109.5 ± 3.4 

12.14 


C 18 CH1 -0) 

109782.2 

5.78 + 0.10 

300.4+ 1.5 

103.0 + 3.5 

52.35 


HNCO(5o,5-4 0 , 4 ) 

109905.8 

0.45 ± 0.79 

300 

100 

4.23 

m, t 

13 CO(l -0) 

110201.4 

28.25 + 0.39 

297.3 ± 1.0 

107.1 ±2.3 

240.09 


CH 3 CN (6a - 5k) 

110383.5 

0.28 ± 0.85 

300 

100 

2.56 

m. t 

C lv O(l -0) 

112359.3 

0.38 ± 0.07 

300 

100 

3.55 

t 

CN(lo,!-Oo,i) 

113191.3 

14.34 + 0.19 

291.5 ± 1.1 

122.2 ± 2.8 

69.10 

m 

CN(1 0 ,2-Oo.i) 

113491.0 

28.63 ± 0.37 

291.5 ± 1.1 

122.2 ± 2.8 

194.80 

m 

NS (3i _i - 2i,i) 

115153.9 

0.73 ±0.13 

300 

100 

6.11 

m, b, t 

12 CO(l -0) 

115271.2 

535.8 + 2.85 

294.3 ± 1.4 

94.4 + 3.1 

3892.23 

b 

Hor38 

115274.4 

- 

300 

100 

9 

b 

NS(3 U -V0 

115556.2 

0.73 + 0.13 

300 

100 

5.55 

m, t 


See caption of Table lA. ll for details. 
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Table A.4: Results for M 51 


Line 

Frequency 

MHz 

I 

K km s^ 1 

Tlsr 
km s 1 

AV 
km s~' 

1 MB 

mK 

Comments 

H 13 CN (1-0) 

86340.2 

0.19 ±0.04 

470 

130 

1.38 

hf, t 

C 2 H(1 -0) 

87316.9 

1.26 ± 0.07 

466.6 ± 5.5 

140.4 ± 12.2 

5.04 

hf 

HCN (1-0) 

88631.8 

5.39 ±0.10 

473.0 ± 1.5 

129.1 ±3.6 

38.75 

hf 

HCO + (1-0) 

89188.6 

2.67 ± 0.07 

473.3 ± 2.2 

132.0 ±5.1 

18.98 


HNC(1 -0) 

90663.6 

1.79 ±0.06 

472.9 ± 2.6 

126.7 ± 6.0 

13.23 


N 2 H + (1 -0) 

93173.7 

0.74 ± 0.04 

466.4 ± 5.5 

146.8 ± 13.0 

4.70 

m 

C 34 S (2-1) 

96412.9 

0.17 ±0.04 

470 

130 

1.24 

t 

CH 3 OH(2 A .-L) 

96741.4 

0.24 ± 0.06 

470 

130 

1.72 

m, t 

CS(2- 1) 

97981.0 

1.04 ±0.05 

464.8 ±4.1 

125.6 ±9.6 

7.74 


SO (2 t - 3 2 ) 

99299.9 

0.25 ± 0.06 

470 

130 

1.80 

t 

C 18 0(1 -0) 

109782.2 

1.61 ±0.09 

474.3 ± 4.7 

133.5 ± 11.1 

11.34 


13 CO(l -0) 

110201.4 

6.68 ± 0.25 

469.7 ± 3.0 

123.9 ±7.1 

50.40 


C lv O(l -0) 

112359.3 

0.30 ± 0.06 

470 

130 

2.14 

t 

CN (l 0 ,i - 0 0>1 ) 

113191.3 

2.55 ± 0.04 

468.4 ± 1.3 

130.0 ±3.4 

12.16 

m 

CN(l 0 , 2 -Oo,i) 

113491.0 

5.12 ±0.07 

468.4 ± 1.3 

130.0 ±3.4 

34.53 

m 

12 CO(l -0) 

115271.2 

50.14 ± 1.53 

468.5 ± 2.5 

120.0 ± 5.7 

379.30 



See caption of Table lA. il for details. 


Table A. 5: Results for NGC 1068 


Line 

Frequency 

MHz 

/ 

Kkms 4 

Tlsr 
km s - 1 

AV 
km s 1 

ry,peak 

1 MB 

mK 

Comments 

o 

1 

Z 

U 

X 

86340.2 

0.78 ± 0.08 

1088.4 ± 11.2 

180 

4.03 

hf 

HCO(l,,o-0,,o) 

86670.8 

0.49 ±0.12 

1100 

240 

1.08 

t, b, hf 

SiO (2 - 1) 

86847.0 

0.60 ± 0.07 

1100 

240 

2.36 

t 

HN I 3 C(1 -0) 

87090.8 

0.19 ±0.05 

1100 

150 

1.20 

t 

C 2 H(1 -0) 

87316.9 

7.80 ±0.10 

1106.4 ± 2.6 

239.8 ±5.3 

18.51 

b, hf 

HNCO(4 0 ,4-3o, 3 ) 

87925.2 

0.57 ± 0.89 

1100 

240 

2.20 

m, t, b 

HCN (1-0) 

88631.8 

23.55 ±0.28 

1112.8 ± 1.6 

215.9 ±3.4 

77.84 

hf 

HCO + (1-0) 

89188.6 

14.36 ±0.17 

1113.7 ± 1.9 

219.3 ±4.2 

56.30 


HOC + (l - 0) 

89487.4 

0.31 ±0.06 

1110 

234.1 ±73.4 

1.23 

t 

HNC(1 - 0) 

90663.6 

7.78 ±0.10 

1112.2 ± 1.9 

225.1 ±4.4 

31.12 


HC,N(10-9) 

90979.0 

0.71 ±0.03 

1091.0 ±5.6 

190 

3.52 


CH 3 CN(5*-4*) 

91987.0 

0.59 ± 0.42 

1100 

240 

2.30 

m, t 

N 2 H + (1 -0) 

93173.7 

1.99 ±0.08 

1100.6 ±5.4 

240 

7.62 

m 

C 34 S (2 - 1) 

96412.9 

0.68 ± 0.08 

1100 

240 

2.64 

t 

C1LOI1 (2/ - L) 

96741.4 

1.55 ±0.09 

1100 

240 

5.93 

m 

CS (2 - 1) 

97981.0 

3.95 ± 0.09 

1092.8 ± 3.3 

212.2 ±7.6 

17.14 

b 

SO(2j-3 2 ) 

99299.9 

0.47 ± 0.06 

1065.1 ± 12.4 

139.7 ±29.0 

3.13 


HC 3 N(11 - 10) 

100076.4 

0.80 ± 0.04 

1091.0 ±5.6 

190 

3.95 


HC 3 N (12 - 11) 

109173.6 

0.85 ± 0.04 

1091.0 ±5.6 

190 

4.19 

t 

C ls O(l -0) 

109782.2 

3.90 ± 0.09 

1100.5 ±4.0 

258.0 ± 9.4 

14.02 


13 CO(l -0) 

110201.4 

13.10 ±0.24 

1100.8 ±2.9 

237.6 ± 5.8 

49.36 


CH 3 CN(6*-5*) 

110383.5 

0.70 ± 0.49 

1100 

240 

2.71 

m, t 

CN(lo,i -Oo.i) 

113191.3 

13.15 ±0.10 

1114.0 ± 1.0 

224.4 ± 2.5 

43.20 

m 

CN(1 0 , 2 -0 0 ,i) 

113491.0 

25.79 ±0.19 

1110.4 ± 1.0 

224.4 ± 2.5 

87.19 

m 

12 CO(l -0) 

115271.2 

81.30 ± 1.01 

1168.1 ± 10.2 

122.5 ± 10.8 

624.39 

m 


CO was fitted with a triple Gaussian, and its parameters refer to the central one. See caption of Table lA. fl for more details. 
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Table A. 6 : Results for NGC 7469 


Line 

Frequency 

MHz 

/ 

Kkms~‘ 

Tlsr 
km s -1 

AV 
km s 1 

~-,peak 

1 MB 

mK 

Comments 

C 2 H(1 -0) 

87316.9 

1.58 ±0.14 

4894.1 ± 17.8 

233.0 ±35.1 

4.19 

hf 

HCN (1 - 0) 

88631.8 

2.36 ±0.13 

4854.3 ± 8.7 

239.2 ± 20.4 

9.18 

hf 

HCO + (1-0) 

89188.6 

2.89 ±0.14 

4847.3 ± 8.0 

247.5 ± 18.7 

10.90 


HNC (1 -0) 

90663.6 

1.08 ± 0.12 

4852.1 ± 18.3 

250 

4.07 


CS (2 - 1) 

97981.0 

0.93 ± 0.09 

4871.0 ± 15.3 

251.2 ±36.0 

3.48 

t 

C ls O(l -0) 

109782.2 

0.42 ± 0.08 

4850 

250 

1.58 

t 

13 CO(l -0) 

110201.4 

2.34 ±0.10 

4849.5 ± 6.9 

245.9 ± 16.1 

8.92 


CN(l 0 ,i - 0 0 j) 

113191.3 

1.92 ±0.07 

4847.3 ± 5.9 

248.0 ± 14.4 

6.33 

m, t 

CN(1 0 , 2 -0 0 ,i) 

113491.0 

3.85 ±0.49 

4847.3 ± 5.9 

248.0 ± 14.4 

14.14 

m 

12 CO(l -0) 

115271.2 

58.12 ±1.51 

4847.0 ±7.6 

236.6 ± 10.1 

208.61 



See caption of Table lA.il for details. 


Table A.7: Results for Arp 220 


Line 

Frequency 

MHz 

/ 

Kkms 4 

Tlsr 
km s ” 1 

AV 

kms " 1 

7-. peak 

1 MB 

mK 

Comments 

H 13 CN (1 - 0) 

86340.2 

1.65 ±0.09 

5364.1 ± 11.1 

328.0 ± 25.6 

4.46 

hf 

H 13 CO + (1 - 0) 

86754.3 

0.40 ± 0.34 

5350 

390 

0.95 

b, t 

SiO (2 - 1) 

86847.0 

1.35 ±0.18 

5350 

363.8 ± 44.0 

3.42 

b 

C 2 H(1 -0) 

87316.9 

3.84 ±0.13 

5317.3 ± 11.5 

388.2 ±32.5 

6.29 

hf 

HCN (1-0) 

88631.8 

11.06 ±0.24 

5339.1 ±6.1 

448.1 ± 13.2 

18.02 

hf 

HCO + (1-0) 

89188.6 

4.62 ±0.21 

5312.8 ± 14.3 

450 

9.13 


HNC(1 - 0) 

90663.6 

7.78 ±0.34 

5326.7 ± 12.7 

378.8 ±27.8 

17.24 


HC,N(10-9) 

90979.0 

3.38 ± 0.09 

5331.4 ±7.4 

386.7 ± 17.4 

8.11 


HC 3 N,v7=l (10_i -9i) 

91202.7 

0.39 ± 1.01 

5331 

390 

1.67 

t 

HC 3 N,v7=l (10 j -9_0 

91333.4 

0.69 ± 1.01 

5331 

390 

1.67 

t 

CH 3 CN(5,-4,) 

91987.0 

1.41 ±0.20 

5350 

386.1 ±83.2 

3.32 

m 

N 2 H + (1 -0) 

93173.7 

2.36 ± 0.42 

5353.0 ± 8.9 

386.3 ± 20.6 

5.42 

m 

C 34 S (2 - 1) 

96412.9 

1.19 ±0.14 

5393.6 ± 35.7 

450 

2.46 

t 

CH 3 OH(2 A ,-l A ) 

96741.4 

0.88 ± 0.09 

5350 

390 

2.09 

m, t 

CS (2 - 1) 

97981.0 

4.15 ±0.07 

5342.6 ± 4.7 

408.1 ± 10.7 

9.14 


HC 3 N(11 - 10) 

100076.4 

4.17 ±0.11 

5331.4 ±7.4 

386.7 ± 17.4 

10.00 


HC 3 N,v7=1 (11, - 10-i ) 

100322.4 

0.99 ± 1.01 

5331 

390 

2.38 


HC 3 N,v 7=1 (11 ! - 10i) 

100466.2 

0.99 ± 1.01 

5331 

390 

2.39 


HC 3 N (12 - 11) 

109173.6 

4.91 ±0.14 

5331.4 ±7.4 

386.7 ± 17.4 

11.78 


HC 3 N,v6=1 (12 - 11) 

109422.0 

1.36 ± 1.01 

5331 

390 

3.28 


HC 3 N.v7=1 (12i — Hi) 

109598.8 

1.37 ± 1.01 

5331 

390 

3.29 


C 18 0(1 -0) 

109782.2 

4.48 ± 0.25 

5330 

450 

9.02 

b 

HNCO(5o,5-4 0 , 4 ) 

109905.8 

1.05 ±0.39 

5344.2 ± 84.4 

367.0 ± 197.3 

2.64 

m, b 

13 CO(l -0) 

110201.4 

4.83 ±0.25 

5330 

450 

9.71 

b 

CH 3 CN(6 A -5 A ) 

110383.5 

1.77 ±0.26 

5350 

386.1 ±83.2 

4.18 

m, b 

CN(lo,i -Oo.i) 

113191.3 

2.98 ± 0.08 

5365.4 ± 10.4 

532.4 ± 19.8 

4.92 

m 

CN(1 0 , 2 -0 0 ,i) 

113491.0 

5.94 ±0.16 

5365.4 ± 10.4 

532.4 ± 19.8 

9.67 

m 

12 CO(l -0) 

115271.2 

94.90 ±0.18 

5336.2 ± 9.3 

429.0 ± 0.9 

207.80 



See caption of Table lA. fl for details. 
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Aladro et al. (2015): A =3 mm line survey of nearby active galaxies 
Table A. 8 : Results for Mrk 231 


Line 

Frequency 

MHz 

I 

Kkms -1 

Vlsr 

kms~' 

AT 
km s 1 

7-, peak 

1 MB 

mK 

Comments 

H 13 CN (1 -0) 

86340.2 

0.46 ±0.10 

12053.5 ± 29.4 

210 

2.02 

hf, t 

C 2 H(1 -0) 

87316.9 

0.57 ±0.11 

12007.0 ±35.1 

210 

1.46 

hf, t 

HCN (1 -0) 

88631.8 

2.06 ±0.12 

11973.2 ± 8.0 

207.8 ± 18.3 

8.86 

hf 

HCO + (1 -0) 

89188.6 

1.34 ± 0.09 

11951.0 ± 10.9 

234.0 ±25.1 

5.23 


HNC(1 - 0) 

90663.6 

0.68 ±0.14 

11975.6 ±28.3 

204.7 ± 65.9 

3.06 


HC,N(10-9) 

90979.0 

0.25 ± 0.06 

12067.4 ± 33.4 

210 

1.10 

t 

CS (2 — 1) 

97981.0 

0.44 ± 0.09 

12048.0 ± 27.4 

210 

1.96 

t 

HC 3 N(11 - 10) 

100076.4 

0.32 ± 0.08 

12067.4 ± 33.4 

210 

1.40 

t 

HC 3 N(12- 11) 

109173.6 

0.38 ± 0.09 

12067.4 ± 33.4 

210 

1.68 

t 

C is O(1 -0) 

109782.2 

0.35 ± 0.06 

12056.6 ± 23.7 

210 

1.57 


l3 CO(l -0) 

110201.4 

0.41 ±0.06 

12086.1 ±21.2 

210 

1.81 


CN(l 0il -0 0 ,i) 

113191.3 

0.31 ±0.21 

12066.7 ± 8.5 

218.7 ± 19.4 

1.32 

m, t 

CN(1 0 .2-Oo,i) 

113491.0 

0.82 ± 0.22 

12066.7 ± 8.5 

218.7 ± 19.4 

3.46 

111 

l2 CO(l -0) 

115271.2 

12.98 ± 1.29 

12038.6 ±7.8 

210 

44.25 



See caption of Table lA. fl for details. 
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